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ABSTRACT 


A user friendly program was developed to calculate fatigue life using Strain- 
Life equations, given either a stress or strain history. Additionally, the material 
parameters and associated stress concentration factors can be varied. Since certain 
material constants, such as cyclic strength coefficients (K’) and strain hardening 
exponents (n’) vary during a material’s fatigue life, the program is capable of either 
keeping them constant or varying them as a function of elapsed cycles. The program 
was then utilized to examine the effects of varying K’ and n’ on the calculated fatigue 
life of aluminum 7075-T6 under a typical flight load history. 
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Le INTRODUCTION 


There are distinctly different approaches used to 
calculate ‘cumulative fatigue damage during the crack 
initiation and crack propagation stages. Since the United 
States Navy, Naval Air Systems Command considers any cracked 
part to be "failed", the focus of this thesis was on crack 
initiation. 

The definition of fatigue damage during the initiation 
phase of a crack’s life is difficult. The damage during the 
crack initiation phase can be related to dislocations, slip 
bands, microcracks, etc, but the phenomena are microscopic and 
are not easily correlated with macroscopic measurements. 
Because of this, the damage surming methods typically used to 
calculate crack initiation are empirical in nature. They 
relate damage to life consumed. Life here refers to the 
physical separation of a small test specimen which is 
subsequently used to approximate crack initiation in larger 
aircraft components, since larger components will have a much 
larger critical crack length. 

The most common method for summing damage is to use the 


linear damage rule, also known as Miner's rule: 


ny 


where n; is the number of cycles applied at a given stress 
level and N; is the total fatigue life for constant amplitude 
loading, at that str-ss level. 

When applying Miner's rule to a variable load history, the 
life used up at each load change has to be calculated. This 
can be done with one of several possible strain-life 


equations; for instance: 


/ 
oO, - 0 
Morrow's: as = —_ (2N,)> + €¢(2N,)¢ (2) 


/ / = 

Of - 9 Of — 5,1? 

Manson-Halford: Aé - 2 aypref =e] (2n,)¢ @) 
Or 


/ 5. 
ad a 
Smith-Watson-Topper: 6 As = ice) 


(4) 
(2N,) *>+04€(2N,) *¢ 





All three of these equations require the local change in 
Strain and either the maximum or mean local stress. These can 
be obtained from relating the far field stress to the local 


stress using Neuber’s empirical rule: 


KiSe = o€ (5) 


Relating stress changes to strain changes can be done with the 


monotonic stress-strain equation: 





or the hysteresis equation: 
2 eg |" (7) 


‘Reference 1) 

Programs are available to solve strain-life equations; 
however, a fatigue life program was needed for research at the 
Naval Postgraduate School, where the user can modify the 
solution algorithms or have access to the prcegram’s source 
code to make changes and explore various facets of the theory. 

For example, to evaluate the effects of varying the cyclic 
strength coefficient and the cyclic strain hardening exponent 
throughout the duration of applied loading, a program named 
FLP was developed. This program solved the three strain-life 
equations, (Equations (2) through (4),) using either fixed or 
varying values for n’ and K’ as specified by the user. 

To develop a realistic sequence of loads from a known 
spectrum, LOADGEN was created. It was used to develop a 
realistic load history for an A-6 based on the three a "g" 
count data developed by LT Rich Walters [Reference 2]. The 
load sequences were processed with the program using various 


strain-life equations. 


II. STRAIN-LIFE COMPUTATION PROGRAM 


A. GENERAL 

The Fatigue Life Program (FLP) was written using 
Microsoft’s QuickBasic. QuickBasic is a relatively simple 
programming language, but an updated and more capable version 
of BASIC, which allows the program to be compiled into an .EXE 
file, executable on any computer using MS-DOS or PC-DOS 
operating system 2.1 or later. 

FLP was designed to be as "user friendly" as possible, 
and to the greatest extent practical, menus were used to set 
the various options. The program was built in a modular 
fashion and documentation was included throughout the program 


to facilitate debugging and later modifications. 


B. MAIN PROGRAM 

The main program controls the general processing flow. It 
starts by establishing certain constants and variables to be 
used throughout the program, and gives initial values to most 
of the user selected options, which could be considered 
default values. Examples of default values would be British 


units and aluminum 7075, which appear as being the units and 


material selected when the program is first started. 


The program then enters a perpetual DO loop wherein it 


sets all the user definable options, reads in the specified 





load sequence, processes the load sequence, and outputs 
processed data to the designated file before returning to the 
user definable options. The loop is exited and the program 
terminated by pressing ESCAPE when the option menus are being 
displayed. 

The main program also sets the video configuration for the 
program, allowing it to automatically select the best mode 
available, and incorporating error trapping sequences to 


indicate when insufficient memory is available. 


C. USER SELECTABLE OPTIONS 
The selection of the various options available to the user 
is accomplished through several menus, each of which is a 
subroutine, or set of subroutines. There are two menus to 
define user selectable options, and a third to input 
properties for new materials. All three of these menus print 
a list of key instructions at the top of the screen, and call 
on another subroutine to update the display screen. 
The first menu to appear is shown in Figure 1 and allows 
the operator to: 
@® choose between using British/U.S. (Brit) units or Standard 
International (SI) units; 
@ choose a local stress concentration factor (K,); 
@® choose a method to calculate the fatigue stress 


concentration factor, either by use of Neuber’s method, 
Peterson’s method, or manual entry; 


@® choose another screen mode if more than one is available; 








@ choose a material from the existing material data base or 
manually enter a new material, which may be saved in the 
material data base; 

@ review the material properties associated with the 
selected material. 

Prior to exiting this menu, the program checks to ensure that 
a value has been either entered or calculated for the fatigue 
stress concentration factor (Kg). If a value for Kg other 
than zero isn’t present, the program will prompt the user for 


one. The program may also be terminated normally from the 


first options menu by pressing ESCAPE. 


sp dey waka deraty Move to next field 

siete ates Move to previous field 
LEFT/RIGHT .... Change field up/down 
Ft .... Display matrial’s parameters 
ENTER .... Start with current values 
ESCAPES cess craleceie cerns siaie Quit Program 


Type of units (SI or British) { Brit ) 


Material { Aluminum 7075-16 as recieved } 
Stress concentration factor (Kt) { 1.000 J 

Method to calculate Kf { Manually Entered ] 

Fatigue stress conc. factor (Kf) { 0.000 J 


Screen Mode ( 12] 





Figure 1 Options Menu 1 


The second menu of user selectable functions, shown in 
Figure 2, appears when the operator exits the first menu by 


pressing RETURN. It allows the operator to: 


@® choose between using stress or strains as input loads: 


@® choose which strain-life equation will be used, either 
Morrow’s, Manson-Halford, or Smith-Watson-Topper; 


@® choose between using either fixed values for the cyclic 
strength coefficient and cyclic strain hardening exponent 
or using values that are a function of elapsed cycles; 

® choose between calculating the cycles to failure, the 
effects of a single load block or processing a 
preprogrammed series of calculations; 


@® select the name of the input file containing the load 
history; 


@® select the name of the output file; 


® activate and deactivate the program’s sound. 


The sound function mentioned above is to assist the 
program’s user when processing a large batch of calculations 
or making a long life calculation. When activated, it will 
beep every time the program finishes processing a load block 
and sound an alarm when completely finished. From the second 
option menu the user has the choice of returning to the first 
option menu by pressing ESCAPE or starting the program’s 


calculation subroutines by pressing RETURN. 


Move to next field 
Move to previous field 
. Change field up/down 
Turn sound to ON 
ENTER Start with current values 
ESCAPE ... Return to previous screen 


Type of inputs (stress or strain) { Stress } 


Equation [ Morrow’s equation } 


Fixed / Varing n’ and K’ { Fixed n' and K’ J] 


Calculation Type { Load blocks to failure J 


Input file name { STRESS.DAT } 


Output file name ( QUTPUT.DAT J 





Figure 2 Option Menu 2 


The other menu used in the program is to facilitate 
operator entry of a new material. Called by the first option 
selection menu, and shown in Figure 3, it displays and allows 
the operator to update the following material properties: 

@ ultimate strength (S,,) 

® yield strength (S,) 
® fatigue yield strength (S,’) 

® strength coefficient (K) 

@ cyclic strength coefficient (K’) 

® strain hardening exponent (n) 

@® cyclic strain hardening exponent (n’) 
@® ductility coefficient (€,) 

® fatigue ductility coefficient (€,') 


@® strength coefficient (d¢) 


@ fatigue strength coefficient (d¢’) 


@® fatigue strength exponent (b) 


@ fatigue ductility exponent (c) 


® endurance strength (S¢) 


@® modulus of elasticity (E) 


Move to next field 
Move to previous field 
LEFT or RIGHT .... Enter a new value 


ENTER ... 
ESCAPE . 


Ultimate Strength 

Yield Strength 

Fatique Yield Strength 
Strength Coeficient 

Cyclic Strength Coefficient 


Strain Hardening Exponent 

Cyclic Strain Hardening Exponent 
Ouctiltity Coefficient 

Fatigue Ductility Coefficient 
Strength Coefficient 

Fatigue Strength Coefficient 
Fatigue strength Exponent 
Fatigue Ductility Exponent 
Endurance Strength 

Modulus of Elasticity 


Return with current values 
... Quit material data entry 


(Su in ksi) 

(Sy in ksi) 

(Sy’ in ksi) 

(K in ksi) 

(K’ in ksi) 

(n) 

(n’) 

(epsilon f) 
(epsilon f’) 
(sigma f in ksi) 
(sigma f’ in ksi) 
(b) 

(c) 

(Sf in ksi) 

(E in ksi) 





Figure 3 Material Entry Menu 


Upon exiting this menu, the user is prompted for the newly 
entered material’s name and gives the option of saving the new 


material in the material data base. 


D. LOAD HISTORY INPUT 

The load history to be processed/analyzed is read into the 
program from the operator designated input file. This is 
accomplished by the subroutine LOADER. Additionally, this 
subroutine determines the number of loads in the sequence and 
whether the first increment is increasing or decreasing. The 
load input file is simply a consecutive list of either 


stresses or strains, one per line, in a DOS compatible format. 


E. CALCULATIONS 

The data processing takes place in the subroutine 
CRUNCHER. In this subroutine the far-field stress or strain 
load sequence is turned into a sequence of far-field stress or 
strain changes, according to what was initially read in. If 


the input file consisted of strain loads, the monotonic 


stress-strain equation, (Equation (6)) and the hysteresis 
equation, (Equation (7),) is used to convert strain changes to 
stress changes. Next using Neuber’s rule, (Equation (5),) the 


far-field stress changes are used to find the local stress 
changes. The local stress changes and either the mean or 
maximum stress are then used in one of the three available 
strain-life equations, (Equations (2) through (4),) to 
determine the number of cycles to failure at the given load 
levels. This life is then used in conjunction with Miner’s 


rule, (Equation (1)) to determine how much damage has 


10 


accumulated, or how much of the fatigue life was consumed, and 
added to a running counter. This process is repeated until 
all the loads are processed, or until failure, which occurs 


when the damage sums to one. 


F. COMPUTATIONAL SUBROUTINES 

Equations(2) through (6) can not be solved directly for 
the needed variables and require an iterative numerical 
solution. Newton's method of successive approximations was 
chosen as the simplest and most efficient method for this 
application. Newton’s method is an iterative method which 
requires the function in question be evaluated at an initial 
estimate. If the estimate produces excessive error, the 
function's derivative is evaluated at the estimated value and 
a new estimate is obtained by subtracting from the old 
estimate the error at the old estimate divided by the 
derivative at the old estimate. [Reference 4] 

Subroutines EQUATIONS1, EQUATIONS2, EQUATIONS3, 
EQUATIONS4, and EQUATIONS5 utilize Newton’s method to solve 
equations (2) through (5) and (7). In EQUATIONS2, EQUATIONS3, 
and EQUATIONS4, which solve the strain-life equations, 
Newton’s method was modified slightly to ensure that variable 
NNf; being solved for was never approximated as being 


negative. (The variable Nnf;, which is raised to a negative 


li 


power, would produce a complex number which the program 
cannot accept.) 

To allow for cycle to cycle variation in K’ and n’ the 
functions xxkf and xxnf were created to redefine the values of 
K' and n’ based on the number of elapsed cycles (calculated by 
the function xxnnfcount). 

Since QuickBasic has no built in function for base 10 
logarithms, the function LOG10 was incorporated to generate 
base ten logarithms from natural logarithms with the 


relationship shown: 


LOg,,(x) = In(io) 


{Reference 5]. 


G. DATA OUTPUT 

Depending on the processing option chosen, the output will 
consist of either a step by step print out of each variable’s 
value for one trip through the input load block, (Figure 4,) 
and a summary of the chosen options, the number of cycles and 
load blocks until failure, or just the output summary, (Figure 
5). The output is appended to the output file so that any 


previous data in the file will not be lost. 


H. CODE VERIFICATION 
To verify the programming code, the FLP was run for 


several short load histories (approximately 10 cycles). The 


first test sequences were constant amplitude loading till 





"Morrow's equation" 

" Fixed n’ and kK’ " 

"Load blocks to failure" 

"input file:", "testi" 

"output file:", "OUTPUT.DAT" 
"blocks:",7 

"i counter:",2482 

"reversal count:",35655 

"life factor:",1.000275410409951 


Figure 4 Sample output summary. 


index Stress deltaStress deltasig i deltaeps NNE 

45 45 45.00 .004369 .15D+09 
20 25 25.00 .002428 .16D409 
90 70 70.00 .006800 .66D405 
60.00 .005826 .14D+06 

20.01 -001942 .18D+09 

40.01 003884 .17D+09 

30.02 .002915 .14D+09 

130.00 .012858 .18D+05 

150.00 -015193 .59D+04 

70.00 .006800 - 200406 

100. .009748 .15D+05 

130. -012858 .58D+04 

70. .006800 .15D+09 

30. .002915 .15D+09 

60. .005826 .12D+09 

50. .004856 .10D+09 

20. .001942 .11D+09 

60. .005826 .11D+09 

117. -011474 .91D404 

44. -004272 .14D+09 

13. .001262 -11D4+09 

132. .013078 -11D+405 

80. .007775 .18D+09 

2. -000875 -17D+09 

Ss. .000487 .20D409 


ecoceoogoccoocoocooe0coococ ece0oc coc c0o0cSo 
og0o0coc00c0c0coecq0ecccoc9o0o0cae0ceo 


Figure 5 Sample step by step output. 


failure and later load sequences were variable amplitude. The 
identical load sequences processed by hand with an HP-48 hand 


held calculator capable of solving complex equations 


numerically. The results were identical to four, usually five 


13 





or six, significant figures. This however didn’t test the 
programs calculation of varying strength coefficient and 
varying strain hardening exponents. To verify that these were 
being correctly calculated the program was run in the 
QuickBasic environment before being compiled into an 
executable, (.EXE) file. When the program is being run in the 
QuickBasic environment, the operator can pause execution at 
any time and, using what’s referred to as the Immediate 
window, execute other instructions [Reference 6]. Using this 
feature and a hand held calculator, it was simple to pause the 


program at random intervals between 1000 and 500,000 elapsed 


cycles to ensure n’ and K’ had correct values. 





III. MATERIAL DATA BASE 


A. CYCLIC PROPERTIES 

In most metals the stress-strain response is altered by 
repeated loading. Depending on the material and it’s initial 
characteristics (i.e., quenched and tempered, or annealed,) a 
metal might either cyclically harden, soften, or have a mixed 
behavior. Typically a material which has a low yield strength 
relative to it’s ultimate tensile strength will undergo strain 
hardening while a hard material cyclically softens. This is 
reflected in the difference between a material’s monotonic 
strength coefficient (K) and its cyclic strength coefficient 
(K’) as well as the strain hardening exponent (n) and the 
cyclic strain hardening exponent (n’). Figure (6) shows the 
effects of cyclic loading on copper. 

Normally, when performing strain-life calculations K’ and 
n’ are used in all the stress-strain calculations. For low to 
medium cycle fatigue, (107 - 10° cycles,) the effect of the 
change in these material properties is undocumented. To 
account for the dynamic nature of these properties, it was 
postulated that they could be expressed as a function of the 
elapsed cycles, allowing their value to be continually updated 


from cycle to cycle. 
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Figure 6 Hysteresis response of copper. (From Ref. 1) 


B. ARCHIVAL DATA 

The program utilizes a separate file MAT.DAT to hold the 
materials available to the program without manual entry by the 
operator. Appendix z shows the contents of the file MAT.DAT. 

The first entry in the material data base indicates the 
number of materials in the original data base. This number 
couldn't be changed without restructuring the data base into 
a series of records so a second file NEWCOUNT was created just 
to record the number of new materials that were added to the 
data base. The number of materials within the data base is 
needed by the program to allow the option menu to cycle 
through all the materials. 

The material properties used in the data base came from 


Metal Fatigue in Engineering by H. O. Fuchs [Reference 3]. 
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The data base contains both the SI and British/American values 


for all properties and when new materials are added to the 
data base both the SI and British/American values are 
recorded. For numerous materials Reference 3 doesn’t include 
values for K, K’ and n. To get around this potential problem, 
when any of these properties are missing the subroutine 
Loadmaterial, which selects the proper material from the data 
base and initializes the variables in the program to the 
appropriate material properties, will estimate the value of 


the missing property using the following relations: 


/ 
or ne n=2 (9) 


(e,) (en) C 





[Reference 1]. 

Review of the selected material’s properties is 
accomplished by pressing Fl while the first option menu is 
displayed. This activates the display used for entering new 
materials, which will display the properties of the material 
that is already loaded. The displayed properties can’t be 


modified to ensure the integrity of the original data base. 
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Iv. LOAD GENERATING PROGRAM LOADGEN 


A. LOAD SPECTRUM CONCEPT 

In order to evaluate the effect of varying n’ and K' on 
strain-life, a realistic load sequence was desired. Readily 
available was the load spectrum data representing the g 
exceedences of a 0.9973 percentile (three sigma) A-6 aircraft, 
Table I, [Reference 2]. 


TABLE 1: THREE o LOAD SPECTRUM 


Date of Corrected 
Oata FLt. Hours 


12/31/89 
1/31/90 
2/28/90 
3131/90 
5/31/90 


7/31/91 
8/31/91 
9/30/91 





The number of exceedences for each «: -yel was entered 
into the program LOADGEN. This program «': © -elect a random 
number and then pick a corresponding 4g, 5g, 6g, or 7g load. 


The random number generator returned a number between 0 and 1, 








which was compared to the percentage of g occurrences greater 


than 7, then greater than 6, and so on, until the random 
number is less than or equal to the percentage of occurrences 
in excess of a certain g level. The algorithm checks to 
ensure that all the g loads at that level haven’t been 
depleted and decrements the counter for the respective bin. 
If all the loads at a particular g level are depleted and a 
random number was generated picking a depleted g level, then 
another random number is called until all the exceedences have 
been picked from the four g levels. 

As each random number was generated, and ag level 
exceedence chosen, the program immediately calls another 
random number to determine the tenths value of the load to be 
inserted into the sequence. For example, if a 5g load was 
picked to be inserted into the sequence, the second random 
number would determine if it were a value of 5.0, 5.1, 5.2, 
etc. up to 5.9. The generated values were then sequenced in 
an output file for temporary storage. With the assumption 
that the load returned to 1g at the end of every cycle, a 


sequence like the one in Figure 7 was generated. 


B. A-6 STRESS SEQUENCE GENERATION 

After the g sequence was developed, it was converted toa 
stress sequence for processing by FLP. Since the exact 
correlation between g load and the stress it produces on the 


A-6 wasn’t known, the standard Navy design criteria was used 
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Figure 7 Cycle by cycle load application 


to calculate the highest acceptable corresponding stress 
values. The Navy's design criteria is such that the yield 
stress can’t be exceeded at the design g limit, and the 
ultimate stress can’t be exceeded at 1.5 times the design g 
limit. Since the design g limit of the A-6 is +6.5, in the 
worst case either a g loading of 6.5 corresponds to the 
material's yield strength, or a g loading of 9.75 corresponds 
to the material’s ultimate strength. Given both the yield and 
ultimate strengths of the material, the program selects the 
limiting case and uses the corresponding g to stress ratio to 
convert the series of g loadings into a series of stress 


loadings. Since the g count information only gives maximums 





and no minimums, the load profile was assumed to return to 1 
g between each of the peaks. 
As the file of g loads was transformed into stress loads, 


the stresses were written sequentially into a file for use by 


FLP. 





V. APPROXIMATION OF K’ AND n’ 


A. EXPERIMENTAL EFFORTS 

To account for the cyclic change in K’ and n’, these 
values were experimentally determined in test specimens whose 
life under a given cyclic strain load had already been 
determined. Specimens made from aluminum 7075-T6, shown in 
Figure 8, were tested in cyclic strain by LT Byron Smith in an 
investigation of the effects of means strain on strain life 
[reference 7]. Under fully reversed strain loading of +.007 
in/in, test specimens were cycled to 10%, 20%, 30%, and 40% of 


their predetermined total life. These specimens, and two 
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Figure 8 Fatigue test specimen 
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specimens without any cyclic preloading, were then subjected 
to a uniaxial stress-strain test to determine the properties 
in question. 

The testing was performed on the Naval Postgraduate 
School, Mechanical Engineering Department’s MTS machine. Data 
correlating to load and displacement were recorded 
continuously throughout the tests. This data was reduced to 
stress and strain, plotted in Figure 9, and then processed to 
determine the stress coefficients and strain hardening 
exponents at the time of testing. 

Two methods of determining the properties in question were 
attempted. In the first the strain hardening exponents (n) 
were calculated using the following relationship from 


reference (1): 
at necking: n= 1n(1 + e) (10) 


where e is the engineering strain. 


Knowing n, K can be determined with the relationship: 


o = K(e,)4 - K=—* (11) 


(e,)” 





where go is the true stress and Ep is the true elastic strain. 


Taking into consideration the following relationships: 


(12) 
€5 = erotal Zs ee 


(13) 
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€-org) = In(l1 + e) (14) 


o = 3(1 +e); (15) 


where ¢€, is the true elastic strain, and S is the engineering 


stress, this expression can be formed: 


S(1 + e) 


(anda +e) - Sesh)’ 


K= 


The resulting values for n’ and K’ appear in Table 2. 


TABLE 2: DATA REDUCTION RESULTS - METHOD 1 


ee 
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The experimentally derived stress coefficient and strain 
hardening data doesn’t lend itself to generating a function 
based on elapsed cycles or used up strain-life. Looking at 
more of LT Smith's work with strain fatigue in aluminum 7075- 
T6é, shows a large deviation in the test results for strain 
life, (Figure 10). Because of the large deviation in the life 
of a specimen, there are large deviations in specified 
percentages of the sample’s life. This indicates the need to 
test a large number of samples to get results for K and n of 


any statistical significance. 
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Figure 9 Stress-Strain Curves 
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Figure 10 Strain-Life Deviation 


B. ESTIMATION OF K’ AND n’ 

Lacking statistically significant test results, for 
purposes of exercising this program both the strength 
coefficient and the strain hardening exponent were postulated 
to be: 

@® the monotonic values when less than 1000 cycles had 
elapsed, 


@® the cyclic values when more than 500,000 cycles had 
elapsed, (half the endurance life,) and 


® a log-normal function of the number of elapsed cycles when 
between 1,000 and 500,000. 
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To calculate these properties between 1000 and 500,000 cycles 
the following equations were used: 


PE eee, iit eee ee 
vari | £09,)(500000) - Log,, (1000) 


(17) 


Joa x) - Log,,(1000)) +n 


K'-K 


Foo. e00000) = tog, aon) toga x) ~ LOg,,{1000)) + K 


For aluminum 7075-T6 these equations are plotted in Figure 11. 
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VI. COMPUTATIONS 


A. COMPARISON PROCEDURE 

To compare the strain-lives calculated by using the cyclic 
properties with those obtained by varying the properties, a 
set of input load files were generated using the three a 


limits of g load data for the A-6 and the LOADGEN program. 


Each file contained approximately 4800 points, representing 


1000 hours of flight time. The subroutine BATCH was setup to 
repeat the strain-life calculation for each of the four load 
files, using each of the three strain-life equations, and the 
four following calculation methods: 

@ fixed n’ and fixed K’ 

@® fixed n’ and variable K’ 

@® variable n’ and fixed K’ 

@® variable n’ and variable K’ 
This amounted to 48 sets of computations, the results of which 


are shown in Table 3. 


B. DISCUSSION OF RESULTS 

The data in Table 3. show a reduction in the calculated 
fatigue life of 15 to 30% when both the strength coefficient 
and the strain hardening exponent were varied from their 


monotonic values to their cyclic values. It is also clear 








chat the strength coefficient has a much greater impact on the 
fatigue life then does the strain hardening exponent. 
Variation of the strength coefficient consistently produced a 
significant drop in the fatigue life, while varying the strain 
hardening exponent very slightly raised the calculated fatigue 
life. 

All three of the strain-life equations used to make the 
calculations behaved the same, though it seemed that the 
respective importance of the strain harding exponent relative 
to the strength coefficient varied between the equations. 

The effect of varying the strength coefficient can be 
anticipated. From Equation(6), a reduced value of kK’ 
translates into greater strain values locally, and greater 
strain means reduced life. However reduced n’ values also 
mean larger strains and consequent reductions in life. 

The effect of varying only the strain hardening could also 
be anticipated, but it’s impact, compared to that of the 
strength coefficient was unexpected. 

It should be noted that these calculations assumed 
constant values for b and c in Equations (2) through (4), 
where in reality they would change. However this example has 
allowed this unique segment of the program to be exercised and 
to gain some preliminary estimates of the effect varying these 


parameters has on fatigue life. 
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Equation 


Morrow's 


Equation 


Smith- 
Watson- 


Topper 


Manson - 


Halford 


TABLE 3: CYCLES TO FAILURE 


Processing 


Option 


eee 


eee eee 
ore 


Input file: 


TESTO | TEST1 | TEST2 | TEST3 


39018 | 38404 | 38543 | 37899 


33933 | 28927 | 32072 








VII. DISCUSSION/CONCLUSIONS 


The FLP originally conceived has been wiitten and appears 
to function well. It possesses flexibity in the manner in 
which it performs fatigue life calculations, using any one of 
three documented strain-life equations which account for mean 
stress effects. The initial exploration done with it supports 
tne concept that fatigue life, especially in the lower end of 
the fatigue spectrum, might be better calculated if changes in 
material properties are accounted for. To do this, however, 
material properties need to be well defined as functions of 
elapsed cycles. The algorithms used in the FLP varied only 
two of the material’s properties as a function of elapsed 
cycles, K’ and n’. Other material properties could be varied 
in addition to the ones currently varied by the FPL. Four 
candidates are: the fracture strength coefficient (a), the 
ductility coefficient (e€,), the fatigue strength exponent (b) 
and the fatigue ductility exponent (c). These coefficients 
have been shown to vary with elapsed cycles as well, and their 
effect should be explored in future tests and incorporated 


into FLP. 
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APPENDIX A FLP PROGRAM LISTING 


A. MAIN PROGRAM 


DEFSNG A-Z 

DEFINT I-P 

DECLARE SUB BATCH () 
DECLARE SUB LOADER () 
DECLARE SUB CRUNCHER (} 
DECLARE SUB EQUATIONS1 
DECLARE SUB EQUATIONS2 
DECLARE SUB EQUATIONS3 
DECLARE SUB EQUATIONS4 
DECLARE SUB EQUATIONSS 
DECLARE SUB OPTMENU1 () 
DECLARE SUB OPTMENU2 () 
DECLARE SUB MATMENU () 
DECLARE SUB UPDATEMENU () 

DECLARE SUB UPDATEMENU2 () 

DECLARE SUB Loadmaterial (index) 

DECLARE SUB NEWMAT (Ky) 

DECLARE SUB OUTPUTER () 

DECLARE SUB HEADER () 

DECLARE SUB DATADUMP () 

DECLARE SUB Petersonkf () 

DECLARE SUB NeuberKf () 

DECLARE SUB GetConfig () 

DECLARE SUB TYPIN (cstring$, valu) 

DECLARE SUB TYPINSTRING (cstring$, stringname$) 

DECLARE SUB Klaxon (Hi%, lows) 

DECLARE FUNCTION xxnf! () 

DECLARE FUNCTION xxKf! () 

DECLARE FUNCTION xxNNfcounté& () 

DECLARE FUNCTION Rotated (Lower, Upper AS INTEGER, present, Inc) 
DECLARE FUNCTION LOG10! (value!) 


wee wee ws 


REM $DYNAMIC 


’ Constants for best avaiiable screen mode 
CONST VGA = 12 

CONST MCGA = 13 

CONST EGA256 = 9 

CONST EGAG64 = 8 

CONST MONO = 
CONST HERC = 3 
CONST CGA = 1 





' User-defined type to hold information about the mode 
TYPE Config 

Sern AS INTEGER 

Colors AS INTEGER 

Atribs AS INTEGER 


XPix AS INTEGER 

YPix AS INTEGER 

TCOL AS INTEGER 

TROW AS INTEGER 
END TYPE 


DIM VC AS Config 


: variables shared throughout the program 


COMMON SHARED YD, i, top, Kt AS SINGLE, Ktf AS SINGLE, option4 AS STRING 
COMMON SHARED matcount, matindex, matname AS STRING, Su, Sy, Syf, K AS SINGLE 
COMMON SHARED Kf AS SINGLE, n AS SINGLE, nf AS SINGLE, equationname AS STRING 
COMMON SHARED epsf, epsff, sigf, sigff, b, c, Sf, SfSu, E, newstuff 

COMMON SHARED flagl AS INTEGER, flag2 AS INTEGER, flag3 AS INTEGER 

COMMON SHARED flag4 AS INTEGER, flag5 AS STRING, batchno AS INTEGER 

COMMON SHARED nkKtype AS STRING, calictype AS STRING, lgNNfcount AS SINGLE 
COMMON SHARED usedlife AS DOUBLE, block AS INTEGER, lasti AS INTEGER 

COMMON SHARED NNfcount AS LONG, xKf AS SINGLE, xnf AS SINGLE 

COMMON SHARED inputfile AS STRING, outputfile AS STRING, sigmax AS SINGLE 
COMMON SHARED laststress AS SINGLE, strainmax AS SINGLE, lasteps AS SINGLE 
COMMON SHARED lastsig AS SINGLE 


size = 5000 


DIM SHARED stress(size) AS SINGLE, strain(size), sig(size) AS SINGLE 
DIM SHARED deltaeps(size) AS SINGLE, NNf(size) AS DOUBLE, sig0 (size) 
DIM SHARED deltaStress(size) AS SINGLE, deltastrain(size) AS SINGLE 
DIM SHARED deltasig(size), eps(size) 


‘ User-defined type to hold information about the selected options (menul) 
TYPE Options1l 


units AS STRING * 4 

material AS INTEGER 

Kt AS SINGLE 

Ktf£ AS INTEGER 

options AS INTEGER 
END TYPE 


DIM Menu AS Optionsl 


' User-defined type to hold information about the selected options (menu2) 
TYPE Options2 

inputs AS STRING * 6 

equation AS INTEGER 


option3 AS INTEGER ‘nKtype selection 
option4 AS INTEGER ‘calculation type 
END TYPE 


DIM menu2 AS Options2 


‘ Error variables to check screen type 
DIM InitRows AS INTEGER, BestMode AS INTEGER, Available AS STRING 
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‘ Initialize variables 


Menu.units = "Brit" 
Menu.material = 21 
Menu.Kt = 1 

7 Menu.Ktf = 1 
Menu.options = 1 


option4 = "Manually Entered " 
. Kt = 1 

Ktf = 0 

menu2.inputs = "Stress" 


menu2.equation = 1 
menu2.option3 = 1 
menu2.option4 = 


equationname = "Morrow’s equation" 
nKtype = " Fixed n’ and kK’ " 
calctype = "Load blocks to failure" 
inputfile = "STRESS.DAT" 

outputfile = "OUTPUT.DAT" 


flag2 = 1 

flag3 = 0 

flag4 = 0 

flagS = "OFF" y initialises sound on 
batchno = 0 

sigmax = 0 


strainmax = 0 

laststress = 0 
lastsig 
lasteps 


ou 
oO 


0 


r Get best configuration and set initial graphics mode to it 
GetConfig 
vC.Secrn = BestMode 


DO ‘ endless do loop for main program (exited from menu 1) 


Selecti n of material/constants 
IF flag4 = 0 THEN 
DO 
CALL OPTMENU1 
CALL OPTMENU2 
LOOP UNTIL flag2 = 1 
END IF 


IF menu2.option4 = 3 THEN 
batchno = batchno + 1 
BATCH 

END IF 


CALL LOADER 
CALL CRUNCHER 


CALL OUTPUTER 





IF flag4 = O THEN 
IF flagS = "OFF" THEN PRINT "Finished Processing..... press any key" 


IF flag5 "OFF" THEN CALL Klaxon(987, 329) 

CLS 
ELSE 

PRINT USING "Finished Processing Batch ##"; batchno J 
END IF 
sigmax = 0 : 


strainmax = 0 
laststress = 0 
lasteps 0 
lastsig 0 


LOOP UNTIL cowscomehome 
END 


‘ Error trap to make program screen independent 


VideoErr: 
SELECT CASE BestMode ' Fall through until something works 
CASE VGA 
BestMode = MCGA 
Available = "12BD" 
CASE MCGA 
BestMode = EGA256 
Available = "12789" 
CASE EGA256 
BestMode = CGA 
Available = "12" 
CASE CG\ 
BestMode = MONO 
Available = "A" 
CASE MONO 
BestMode = HERC 
Available = "3" 
CASE ELSE 
PRINT "Sorry. Graphics not available. " 
END 
END SELECT 
RESUME 


‘ Trap to detect 64K EGA 


EGAErr: 
BestMode = EGA64 
Available = "12789" 


RESUME NEXT 


' Trap to detect insufficient memory 
MemErr: 

LOCATE 22, 1 

PRINT "Out of memory" 

RESUME NEXT 





‘ Trap to determine initial number of rows so they can be restored 
RowErr: 


IF InitRows = 50 THEN 
InitRows = 43 
RESUME 

ELSE 


InitRows = 25 
RESUME NEXT 
END IF 





B. SUBROUTINE BATCH 


SSESSsrstSs SSS SSaesSesrSSeeSrse= BATCH west sesesssseesesesscoeesestqeSsesoqqeqe=ersF= 
: Subroutine to access a file for repeated runs 
' (currently set to process 4 file, with 3 eqns, 4 methods 


SUB BATCH 
SHARED menu2 AS Options2 


IF batchno = 1 THEN 
inputfile = "testaa" 
nKtype = " Fixed n’ and K’ " 
menu2.option3 = 1 
menu2.equation = i 
equationname = "Morrow’s Equation" 
ELSEIF batchno = 2 THEN 
inputfile = "testaa" 
nKtype = "Variable n’ and K’" 
menu2.option3 = 2 
menu2.equation = 1 


equationname = "Morrow’s Equation" 
ELSEIF batchno = 3 THEN 

inputfile = "testaa" 

nKtype = "Variable n’ and fixed K’ " 


menu2.option3 = 3 

menu2.equation = 1 

equationname = "Morrow’s Equation" 
ELSEIF batchno = 4 THEN 

inputfile = "testaa" 

nKtype = "Fixed n’ and variable K’" 

menu2.option3 = 4 

menu2.equation = 1 


equationname = "Morrow’s Equation" 
ELSEIF batchno = 5 THEN 

inputfile = "testbb" 

nKtype = " Fixed n’ and kK’ " 


menu2.option3 = 1 
menu2.equation = 1 
equationname = "Morrow’s Equation" 
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ELSEIF batchno = 6 THEN 
inputfile = "testbb" 
nKtype = "Variable n’ and K’" 
menu2.option3 = 2 
menu2.equation = 1 


equationname = "Morrow’s Equation" 
ELSEIF batchno = 7 THEN 

inputfile = "testbb" 

nKtype = "Variable n’ and fixed K’ " 


menu2.option3 = 3 

menu2.equation = 1 

equationname = "Morrow's Equation" 
ELSEIF batchno = 8 THEN 

inputfile = "testbb" 

nKtype = "Fixed n’ and variable K‘'" 

menu2.option3 = 4 

menu2.equation = 1 

equationname = "Morrow's Equation" 
ELSEIF batchno = 9 THEN 

inputfile = "testec" 

nKtype = " Fixed on’ and K’ " 

menu2.option3 = 1 

menu2.equation = 1 


equationname = "Morrow’s Equation" 
ELSEIF batchno = 10 THEN 

inputfile = "testca" 

nKtype = "Variable n’ and K'" 


menu2.option3 = 2 
menu2.equation = 1 


equationname = "Morrow’s Equation” 
ELSEIF batchno = 11 THEN 
inputfile = "testcc" 


nKtype = "Variable n’ and fixed K’ " 

menu2.option3 = 3 

menu2.equation = 1 

equationname = "Morrow’s Equation" 
ELSEIF batchno = 12 THEN 

inputfile = "testcc" 

nKtype = "Fixed n’ and variable K’'" 

menu2.option3 = 4 

menu2.equation = 1 

equationname = "Morrow’s Equation" 
ELSEIF batchno = 13 THEN 

inputfile = "testdd" 

nKtype = " Fixed n’ and K’ " 

menu2.option3 = 1 

menu2.equation = 1 


equationname = "Morrow’s Equation" 
ELSEIF batchno = 14 THEN 
inputfile = "testdd" 


nKtype = "Variable n’ and K’" 
menu2.option3 = 2 
menu2.equation = 1 


equationname = "Morrow’s Equation" 
ELSEIF batchno = 15 THEN 
inputfile = "testdd" 


nKtype = "Variable n’ and fixed K’ " 
menu2.option3 = 3 

menu2.equation = 1 

equationname = "Morrow’s Equation" 
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ELSEIF batchno = 16 THEN 

inputfile = "testdd" 

nKtype = "Fixed n’ and variable K’" 

menu2.option3 = 4 

menu2.e@quation = 1 

equationname = "Morrow’s Equation" 
ELSEIF batchno = 17 THEN 

inputfile = "testaa" 

nKtype = " Fixed n’ and K’ " 

menu2.option3 = 1 

menu2.equation = 2? 

equationname = "Smith-Watson-Topper" 
ELSEIF batchno = 18 THEN 

inputfile = "testaa" 

nKtype = "Variable n’ and K’" 

menu2.option3 = 2 

menu2.equation = 2 


equationname = "Smith-Watson-Topper" 
ELSEIF batchno = 19 THEN 
inputfile = "testaa” 


nKtype = "Variable n’ and fixed K’ " 
menu2.option3 = 3 
menu2.equation = 2 


equationname = "Smith-Watson-Topper" 
ELSEIF batchno = 20 THEN 

inputfile = "testaa" 

nKtype = "Fixed n’ and variable K’" 


menu2.option3 = 4 

menu2.equation = 2 

equationname = "Smith-Watson-Topper" 
ELSEIF batchno = 21 THEN 

inputfile = "testbb" 

nKtype = " Fixed n’ and K’ " 

menu2.option3 = 1 

menu2.equation = 2 


equationname = "Smith-Watson-Topper" 
ELSEIF batchno = 22 THEN 
inputfile = "testbb" 


nKtype = "Variable n’ and K’" 
menu2.option3 = 2 
menu2.equation = 2 


equationname = "Smith-Watson-Topper” 
ELSEIF batchno = 23 THEN 
inputfile = "testbb" 


nKtype = "Variable n’ and fixed K’ " 
menu2.option3 = 3 
menu2.equation = 2 


equationname = "Smith-Watson-Topper" 
ELSEIF batchno = 24 THEN 

inputfile = "testbb" 

nKtype = "Fixed n’ and variable K’" 


menu2.option3 = 4 

menu2.equation = 2 

equationname = "Smith-Watson-Topper" 
ELSEIF batchno = 25 THEN 

inputfile = "testcc" 

nKtype = " Fixed n’ and K’ " 

menu2.option3 = 1 

menu2.equation = 2 

equationname = "Smith-Watson-Topper" 
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ELSEIF batchno = 26 THEN 

inputfile = "testcec" 

nKtype = "Variable n‘ and K‘" 

menu2.option3 = 2 

menu2.equation = 2 

equationname = "Smith-Watson-Topper" . 
ELSEIF batchno = 27 THEN 

inputfile = "testcc” 

nKtype = "Variable n’ and fixed K’ " 

menu2.option3 = 3 

menu2.equation = 2 


equationname = "Smith~-Watson-Topper" 
ELSEIF batchno = 28 THEN 

inputfile = "testcc" 

nKtype = "Fixed n’ and variable K’" 


menu2.option3 = 4 

menu2.equation = 2 

equationname = "Smith-Watson-Topper" 
ELSEIF batchno = 29 THEN 

inputfile = "testdd" 

nKtype = " Fixed n’ and K’ " 

menu2.option3 = 2 

menu2.equation = 1 


equationname = "Smith-~Watson-Topper"” 
ELSEIF batchno = 30 THEN 
inputfile = "testdd" 


nKtype = "Variable n’ and K’" 

menu2.option3 = 2 

menu2.equation = 2 

equationname = "Smith-Watson-Topper" 
ELSEIF batchno = 31 THEN 

inputfile = “testdd" 

nKtype = "Variable n’ and fixed K’ " 

menu2.option3 = 3 

menu2.equation = 2 

equationname = "Smith-Watson-Topper" 
ELSEIF batchno = 32 THEN 

inputfile = "testdd" 

nKtype = "Fixed n’ and variable K’”" 

menu2.option3 = 4 

menu2.equation = 2 

equationname = "Smith-Watson-Topper” 
ELSEIF batchno = 33 THEN 

inputfile = "testaa" 

nKtype = " Fixed n’ and K’ " 

menu2.option3 = 1 

menu2.equation = 3 


equationname = "Manson-Halfora" 
ELSEIF batchno = 34 THEN 
inputfile = "testaa" 


nKtype = "Variable n’ and K’" 
menu2.option3 = 2 
menu2.equation = 3 


equationname = "Manson-Halford" 
ELSEIF batchno = 35 THEN 
inputfile = "testaa" 


nKtype = "Variable n’ and fixed K’ " 
menu2.option3 = 3 

menu2.equation = 3 

equationname = "Manson-Halford" 





ELSEIF batchno = 36 THEN 

inputfile = "testaa" 

nKtype = "Fixed n’ and variable 

menu2.option3 = 4 

menu2.equation = 3 

equationname = "Manson-Halford" 
ELSEIF batchno = 37 THEN 

inputfile = "testbb" 

nKtype = " Fixed n’ and K’ " 

menu2.option3 = 1 

menu2.equation = 3 

®quationname - "Manson-Halford" 
ELSEIF batchno = 38 THEN 

inputfile = "testbb" 

nKtype = "Variable n’ and k’" 

menu2.option3 = 2 

menu2.equation = 3 

equationname = "Manson-Halfora" 
ELSEIF batchno = 39 THEN 

inputfile = "testbb" 

nKtype = "Variable n’ and fixed 

menu2.option3 = 3 

menu2.equation = 3 

equationname = "Manson-Halford" 
ELSEIF batchno = 40 THEN 

inputfile = "testbb" 

nKtype = "Fixed n’ and variable 

menu2.option3 = 4 

menu2.equation = 3 

equationname = "Manson-Halford" 
ELSEIF batchno = 41 THEN 

inputfile = "testcec" 

nKtype = " Fixed n’ and K’ " 

menu2.option3 = 1 

menu2.equation = 3 


equationname = "Manson-Halford" 
ELSEIF batchno = 42 THEN 
inputfile = "testcc” 


nKtype = "Variable n’ and K’" 
menu2.option3 = 2 
menu2.equation = 3 


equationname = "Manson-Halford" 
ELSEIF batchno = 43 THEN 
inputfile = "testcc" 


Ke " 


K’ 


K’ w 


nKtype = "Variable n’ and fixed XK’ 


menu2.option3 = 3 

menu2.equation = 3 

equationname = "Manson-Halford" 
ELSEIF batchno = 44 THEN 

inputfile = "testcc" 


nKtype = "Fixed n’ and variable K’" 


menu2.option3 = 4 

menu2.equation = 3 

equationname = "Manson-Halford" 
ELSEIF batchno = 45 THEN 

inputfile = "testdd" 

nKtype = " Fixed n’ and K’ " 

menu2.option3 = 1 

menu2.equation = 3 

equationname = "Manson-Halford" 
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ELSEIF batchno = 46 THEN 
inputfile = "testdd" 
nKtype = "Variable n’ and K’" 
menu2.option3 = 2 
menu2.equation = 3 


equationname = "Manson-Halford" : 
ELSEIF batchno = 47 THEN 
inputfile = "testdd" 


nKtype = "Variable n’ and fixed K’ " 

menu2.option3 = 3 

menu2.equation = 3 

equationname = "Manson-Halford" 
ELSEIF batchno = 48 THEN 

inputfile = "testdd" 

nktype = "Fixed n’ and variable K’" 

menu2.option3 = 4 

menu2.equation = 3 


equationname = "Manson-Halford" 
flag4 = 0 

END IF 

END SUB 


C. SUBROUTINE CRUNCHER 


‘  gubroutine to calculate Strian-Life using the various methods 


SUB CRUNCHER 


DIM lastload AS SINGLE 
SHARED menu2 AS Options2 


OPEN outputfile FOR APPEND AS #7 


‘initialize variables 


usedlife = 0 
block = -1 
lastload = 0 
DO ’ loop until life is used up if single block option not chosen 
block = block + 1 
IF block = 0 THEN ’calculate delta loads and print output header if 
appropriate 
IF menu2.option4 = 2 THEN HEADER ‘’ output headers 
IF menu2.inputs = "Stress" THEN 


calculation of farfield deltastress array 
FOR ii = 1 TO top 
deltaStress(ii) = ABS(stress(ii) - lastload) 
lastload = stress (ii) 
NEXT ii 
ELSE 
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, calculation of far field deltastrain array 
FOR ii = 1 TO top 


deltastrain(ii) = ABS(strain(ii) - lastload) 
lastload = strain(ii) 
NEXT ii 
END IF 
END IF 


’ 


FOR i = 1 TO top 


load blocks processed 


, set properties based on elapsed cycles 
NNfcount = xxNNfcount 
lgNNfcount = LOG10(CSNG(NNfcount) ) 
xKf = xxKf 
xnf = xxnf 
: turn farfield strains into stresses 
IF menu2.inputs = "Strain" THEN CALL EQUATIONSS 


£ find local Stress associated with farfield stress 
CALL EQUATIONS1 
: Calcalate mean local stress 
sig(i) = lastsig + ({-1) * (flagl + i)) * deltasig(i) 
sigO(i) = (lastsig + sig(i)) / 2 
lastsig = sig(i) 
' finding epsilon and delta epsilon 
IF sig(i) > sigmax THEN 
f mopnotonic equation to find local strains 
eps(i) = sig(i) / E + (sig(i) / (xKf)) * (1 / xnf) 
deltaeps(i) = ABS(eps(i) - lasteps) 
lasteps = eps(i) 
sigmax = sig(i) 
ELSE 
y hysteresis equation to find local strains 
deltaeps(i) = deltasig (i) /E+2 * (deltasig(i) / (2 * xKf)) 


(1 / xnf) 
eps(i) = lasteps + ((-1) * (flagl + i)) * deltaeps (i) 
lasteps = eps(i) 
END IF 


. pick between various strain-life eqns 
SELECT CASE menu2.equation 


CASE 1 ’ Morrow’s equation 
CALL EQUATIONS2 

CASE 2 ‘’ Smith Watson Topper 
CALL EQUATIONS3 

CASE 3 *’ Manson Halford 
CALL EQUATIONS4 

CASE ELSE 

END SELECT 


. print running output to file 
IF menu2.option4 = 2 THEN DATADUMP 





IF usedlife > 1 THEN 
lasti =i 
EXIT FOR 
END IF 
NEXT i 


IF flagS = "OFF" THEN BEEP 
LOOP UNTIL usedlife >= 1 OR (menu2.option4 = 2 AND block = i) 


CLOSE #7 


END SUB 


D. SUBROUTINE DATADUMP 


SUB DATADUMP 


SHARED menu2 AS Options2 


IF menu2.inputs = "Stress" THEN 


PRINT #7, USING " ### HEH ### ### #2 H#HH HH 
HHHH HH OF HHHHHH HL HHRHHH HHO AAA OG; 
stress(i); deltaStress({i); deltasig(i); sig(i); 
sigO(i); deltaeps(i); eps(i); NNf (i) 


ELSE 
PRINT #7, USING " ### #.### 8.844 #44 .384 HHH HH RHRH HH 
HEHE HH RLHHHHHH FL HHHHHH 4. HECT**" OA; 
strain(i); deltastrain(i); deltaStress (i) ; 
deltasig(i); sig({i); sigQ(i); deltaeps (i); 
eps(i); NNf£ (i) 





E. SUBROUTINE EQUATIONS1 


subroutine to find delta sigma using Newton’s method 
(based on Neuber’s Relation) 


SUB EQUATIONS1 

DIM Yprime AS DOUBLE 
Jeltasig(i) = 1 
Loopcount = 0 


DO 
loopcount = loopcount + 1 
Y = deltasig(i) * (deltasig(i) / (2 * E) + (deltasig(i) / (2 * xKf)) 
* (1 / xnf)) - (Ktf£ * 2) * deltaStress(i) * (deltaStress(i) / (2 * E) 
+ (deltaStress(i) / (2 * xKf)) * (1 / xnf)) 
IF ABS(Y) > .0000001 THEN 
Yprime = deltasig(i) / E + (({xnf + 1) / xnf) * (deltasig(i) / (2 * xKf)) 
“(1 / xnf) 
deltasig(i) = deltasig(i) - Y / Yprime 
END IF 


LOOP UNTIL (ABS(Y) <= .0000001) OR (loopcount = 10000) 


END SUB 


F. SUBROUTINE EQUATIONS2 


‘ seecsseersrsseSeeeSeeeeeeze=e= EQUATIONS2 th ee 7 ee ee ye ee 
‘ subroutine to evaluate Morrow’s Strain-Life equation 

SUB EQUATIONS2 

REM SDYNAMIC 


DIM Yprime AS DOUBLE, Y AS DOUBLE, loopcount AS LONG 


NN£(i) = 10000 
loopcount = 0 


DO 
loopcount = loopcount + 1 
Y = -deltaeps(i) / 2 + ((sigff - sigO(i)) / E) * (2 * NN£(i)) * b + epsff * 
(2 * NNE(i)) * c 


IF ABS(Y) > .0000000001# THEN 
Yprime = (b * (sigff - sigO(i)) / BE) * (2 * b) * (NNE£(i)) * (b- 1) + 
* epsff * (2 * c) * (NNE(i)) * (ce - 1) 
IF (Y / Yprime < NNf£(i)) THEN 


NNf(i) = NN£f(i) - Y / Yprime 
ELSE 

NN£(i) = NNE(i) / 2 
END IF 


END IF 








LOOP UNTIL (ABS(Y) <=.0000000001#) OR (NN£(i) > 100000000) OR (loopcount=10000) 


usedlife = usedlife + 1 / NNF£ (i) 
IF loopcount = 10000 THEN PRINT "*" ‘indication of convergence dificulties 


END SUB 


G. SUBROUTINE EQUATIONS3 


‘ subroutine to evaluate Smith Watson Topper Strain-Life equation 


SUB EQUATIONS3 
REM SDYNAMIC 


DIM Yprime AS DOUBLE 


NN£(i) = 100 
loopcount = 0 


DO 
loopcount = loopcount + 1 
Y = -(sigQ0(i) + deltasig(i) / 2) * (deltaeps(i) / 2) + ((sigff * 2) / E) * (2 
* NN£(i)) * (2 * b) + sigff * epsf£ * (2 * NN£(i)) * (b +c) 


IF ABS(Y) > .0000000001# THEN 
Yprime = (2 * b * (sigff ~ 2) / E) * (2 * (2 * b)) * NNE£(i) * (2 * b - 1) 
+ (b + c) * sigff * epsff * (2 * (b + c)) * NN£(1) * (b+ c¢ - 1) 
IF (Y / Yprime < NNf£(i)) THEN 
NN£(i) = NN£(i) - Y / Yprime * .5 
ELSE 
NNf (i) = NN£(i) / 2.5 
END IF 
END IF 
LOOP UNTIL (ABS(Y) <= .0000000001#) OR (NN£(i) > 100000000) OR (loopcount = 
10000) 


usedlife = usedlife + 1 / NN£ (i) 
IF loopcount = 10000 THEN PRINT "*" ‘indication of convergence difficulties 


END SUB 
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H. SUBROUTINE EQUATIONS4 


' subroutine to evaluate Manson-Halford’s Strain-Life equation 


SUB EQUATIONS4 
REM S$DYNAMIC 


DIM Yprime AS DOUBLE 


NN£(i) = 100 
loopcount = 0 


DO 
loopcount = loopcount + 1 
Y = -deltaeps(i) / 2 + ((sigff - sigO(i)) / E) * (2 * NN£(i)) * b + epsff * 
((sigff - sigO(i)) / sigff) * (b / c) * (2 * NNE(i)) * ¢ 
IF ABS(Y) > .0000000001# THEN 
Yprime = (b * (sigff - sigO0(i)) / E) * (2 * b) * NN£(i) * (b - 1) +c * 
epsff * ((sigff - sigO(i)) / sigff) * (b/c) * (2 * c) * NNE(i) 


“(ce - 1) 
IF (Y / Yprime < NNf(i)) THEN 
NN£ (i) = NN£(i) - Y / Yprime * .5 
ELSE 
NNf (i) = NN£(i) / 2.5 
END IF 
END IF 


LOOP UNTIL (ABS(Y) <=.0000000001#) OR (NN£(i) > 100000000) OR (loopcount=10000) 


usedlife = usedlife + 1 / NNf (i) 


IF loopcount = 10000 THEN PRINT "*" ‘indication of convergence dificulties 


END SUB 


I. SUBROUTINE EQUATIONSS5 


‘ gubroutine to find the far-field stress associated with the 
. far-field strain inputs 


SUB EQUATIONSS 
DIM Yprime AS DOUBLE 


deltaStress(i) = 0 
loopcount = 0 


IF strain(i) < strainmax THEN 
DO 
loopcount = loopcount + 1 
Y = -deltastrain(i) + deltaStress(i) / E + 2 * (deltaStress / (2 * xKf)) 


“(1 / »nf) 
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IF ABS{Y) > .0001 THEN 
Yprime = 1/ E+ (1 / (2 * xKf)) * (1 / xnf) * (2 / xnf) * deltaStress 
“(1 / xnf - 1) 
deltaStress(i) = deltaStress(i) - Y / Yprime 
END IF 
LOOP UNTIL (ABS(Y) <= .0001) OR (loopcount = 1000) 


' calcalate far field stress from change in stress 
stress(i) = laststress + ((-1) * (flagl + i)) * deltaStress (2) 
laststress = stress (i) 
ELSE 
DO 
loopcount = loopcount + 1 
Y = -strain(i) + stress(i) / E + (stress / xKf) ~ (1 / xnf) 
IF ABS(Y) > .0001 THEN 
Yprime = 1 / E + (1 / (xKf)) * (1 / »mf) * (1 / »xmf) * stress * (1 
/ xnf - 1) 
stress(i) = stress(i) - Y / Yprime 
END IF 
LOOP UNTIL (ABS(Y) <= .0001) OR (loopcount = 1000) 


‘ cal. farfield stress change and reset strainmax 
strainmax = strain(i) 


deltaStress(i) = ABS(stress({i) - laststress) 
laststress = stress (i) 
END IF 
END SUB 


Jd. 


‘ 


SUBROUTINE GetConfig 


SSBPSSessesSsSS2Seeseeersesereeeq GetConfig SPSRSsesSSSSSSSSSESeesqsssesesszuz 
Get the starting number of lines and the video adapter. 


SUB GetConfig STATIC 
SHARED InitRows AS INTEGER, BestMode AS INTEGER, Available AS STRING 


’ Assume 50 line display and fall through error 
’ until we yet the actual number 

InitRows = 50 

ON ERROR GOTO RowErr 

LOCATE InitRows, 1 


’ Assume best possible screen mode 
BestMode = VGA 
Available = "12789BCD" 


ON ERROR GOTO VideoErr 

‘ Fall through error trap until a mode works 

SCREEN BestMode 

' If EGA, then check pages to see whether more than 64K 
ON ERROR GOTO EGAErr 

IF BestMode = EGA256 THEN SCREEN 8, , 1 
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ON ERROR GOTO 0 


’ 


Reset text mode 
SCREEN 0, , 0 
WIDTH 80, 25 


END SUB 


K. SUBROUTINE HEADER 


eee ee ees HEADER SrrSssereSsrrtstssesSSSsssSssserrenesrzrssz 


SUB HEADER 
SHARED menu2 AS Options2 


IF menu2.inputs = "Stress" THEN 
PRINT #7, "index Stress @GeltaStress deltasig 


deltaeps eps NNf " 
ELSE 
PRINT #7, "index strain dstrn deltaStress deltasig 
deltaeps eps NNf " 
END IF 
END SUB 


L. SUBROUTINE Klaxon 


sig 


sig 


sigOd 


sigd 


SrSeSsrsessesseseseerrsesssessere=sz Klaxon SSSseesrsesss2Sseq2eseesressessseeszerzeqsz 


SUB Klaxon (Hit, lows) STATIC 
DO WHILE INKEYS = "" 
SOUND Hi, 5 
SOUND low, 5 
LOOP 


END SUB 
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M. SUBROUTINE LOADER 


REM $DYNAMIC 


" subroutine to read in the array of loads from the file "STRESS.DAT" 
: or "STRAIN.DAT" as appropriate 

m also -- "top" is set to the highest index in the load array and 
-- "flagl" is set to indicate if the loads are initially 

4 increasing or decreasing 


SUB LOADER 
SHARED menu2 AS Options2 


OPEN inputfile FOR INPUT AS #10 
j= 


IF menu2.inputs = "Stress" THEN 


DO UNTIL EOF (10) 
INPUT #10, stress(j) 
j=j+l 
LOOP 


t set "top" to the number of loads: 
top =j-l 


, set "flagi": (indicating first cycle 
IF stress(1) > 0 THEN ‘ in tension or compression) 
flagi = 1 
ELSE 
flagi = 0 
END IF 
ELSE 
DO UNTIL EOF (10) 
INPUT #10, strain(3) 
3 = 3j + 1 
LOOP 
. set "top” to the number of loads: 
top =j- 1 
set "flagl": (same as above) 
IF strain(1) > 0 THEN 
flagl = 1 
ELSE 
flagi = 0 
END IF 
END IF 
CLOSE #10 
END SUB 
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N. SUBROUTINE Loadmaterial 


REM $STATIC 


=SSSscSS Sess SS Sse sSSsesees = Loadmaterial esSSs5SSSSSSsSS2SSsSSS2==S=>55= 


: finds the selected material in the data base and loads the 
4 Material values for the selected material 
af also calculates any critical missing values it can 


SUB Loadmaterial (index) 

‘OPEN "A:\MAT.DAT" FOR INPUT AS #3 
“OPEN "B:\MAT.DAT" FOR INPUT AS #3 
OPEN "MAT.DAT" FOR INPUT AS #3 
INPUT #3, count 

DO UNTIL index = matindex 


INPUT #3, matindex, matname$, Su, Sy, Syf, K, Kf, n, nf, epsf, epsff, 
sigff, b, c, Sf, SfSu, E 


LOOP 
CLOSE #3 

IF n = 0 AND b <> 0 AND c <> 0 THENn = b/c 

IF nf = 0 AND b <> 0 AND c <> O THEN nf = b/e 

IF K = 0 AND n <> 0 AND epsf <> 0 THEN K = sigf / (epsf * n) 

IF Kf = 0 AND nf <> 0 AND epsff <> 0 THEN Kf = sigff / (epsff * nf) 


END SUB 


O. FUNCTION LOG10 


FUNCTION LOG10! (value) 
LOGi0 = LOG(value) / LOG(10) 


END FUNCTION 
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sigf, 


P. SUBROUTINE MATMENU 


* seaeseseretesesens=aeseteaesS==== MATMENU iF — a 
, subroutine to update the screen when entering or viewing materials 
in the materials database 


CONST ENTER = 13, ESCAPE = 27 

CONST DOWNARROW = 80, UPARROW = 72, LEFTARROW = 75, RIGHTARROW = 77 
CONST COL1 = 10, COL2 = 50, ROW = 7 

CONST COL3 = 7, COL4 = 62 

CONST Fields = 14 


SHARED Menu AS Optionsl 
DIM forceunit AS STRING, Fld AS INTEGER 


CLS 

IF flag3 = 0 THEN 
‘ Display key instructions 
LOCATE 1, COL1 


PRINT "UP! cis d Sate dceigasiete Move to next field" 
LOCATE 2, COL1 

PRINT "DOWN ........ Move to previous field" 
LOCATE 3, COL1 

PRINT "LEFT or RIGHT .... Enter a new value" 
LOCATE 4, COL1 

PRINT "ENTER ... Return with current values" 
LOCATE 5, COL1 

PRINT "ESCAPE .... Quit material data entry" 

ELSE ‘display material’s name at top of screen 

LOCATE 2, COL1 

PRINT "ESCAPE ....to return to previous previous menu" 


LOCATE 5, COL1 
PRINT USING "Data base parameters for: &"; matname 


END IF 
IF (Menu.material <= 51) THEN 

forceunit = "ksi)" 

ELSE 

forceunit = "MPa)" 

END IF 
' Display fields 
LOCATE ROW, COL3: PRINT "Ultimate Strength (Su in "; 
forceunit 


LOCATE ROW, COL4: PRINT USING "( #### ]"; Su; 


LOCATE ROW + 1, COL3: PRINT "Yield Strength (Sy in "; 
forceunit; 
LOCATE ROW + 1, COL4: PRINT USING "[ #### ]"; Sy; 


LOCATE ROW + 2, COL3: PRINT "Fatigue Yield Strength (Sy’ in "; 
forceunit; 


LOCATE ROW + 2, COL4: PRINT USING "[ #### ]"; Syf; 





LOCATE 
LOCATE 
LOCATE 
LOCATE 


LOCATE 
LOCATE 


LOCATE 
LOCATE 


LOCATE 
LOCATE 


LOCATE 
LOCATE 


LOCATE 
LOCATE 
LOCATE 
LOCATE 


LOCATE 
LOCATE 


LOCATE 
LOCATE 


LOCATE 
LOCATE 
LOCATE 


LOCATE 


END SUB 


ROW 
ROW 
ROW 
ROW 


ROW 
ROW 


ROW 
ROW 


ROW 
ROW 


ROW 
ROW 


ROW 
ROW 
ROW 
ROW 


ROW 
ROW 


ROW 
ROW 


ROW 
ROW 
ROW 


ROW 


+ 


+ 
+ 


3, COL3: PRINT "Strength Coeficient (K in "; 
forceunit; 

3, COL4: PRINT USING "[ #### ]"; K; 

4, COL3: PRINT "Cyclic Strength Coefficient (K’ in "; 
forceunit; 

4, COL4: PRINT USING "[( #### ]"; Kf; 

5S, COL3: PRINT "Strain Hardening Exponent (n)"; 

S, COL4: PRINT USING "{ #.## ]"; n; 

6, COL3: PRINT "Cyclic Strain Hardening Exponent (n’)"; 

6, COL4: PRINT USING "[{ #.## ]"; nf; 

7, COL3: PRINT "Ductility Coefficient (epsilon f)"; 

7, COL4: PRINT USING "[ #.## J"; epsf; 


+ 8, COL3: PRINT "Fatigue Ductility Coefficient 


+ 


+ 


+ 


+ 


+ 


+ 
+ 


8, COL4: PRINT USING "[ #.## ]"; epsff; 

9, COL3: PRINT "Strength Coefficient 

9, COL4: PRINT USING "[ #### ]"; sigf; 

10, COL3: PRINT "Fatigue Strength Coefficient 
10, COL4: PRINT USING "[{ #### ]"; sigff; 


11, COL3: PRINT "Fatigue strength Exponent 
11, COL4 - 1: PRINT USING "[ ##.## ]"; b; 


12, COL3: PRINT "Fatigue Ductility Exponent 
12, COL4 - 1: PRINT USING "[ ##.## ]"; c; 


13, COL3: PRINT "Endurance Strength 
13, COL4 - 1: PRINT USING "[ ##### ]"; Sf; 
14, COL3: PRINT "Modulus of Elasticity 


14, COL4 - 1: PRINT USING "[ ##### ]"; E; 
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(epsilon f')"; 


(sigma f in "; 
forceunit 


(sigma f’ in "; 
forceunit 


(b) " : 


(c) Mos 


(Sf in "; 
forceunit 


(E in "; 
forceunit 


Q. SUBROUTINE NeuberKf 


a ee ee ee eee NeuberKf SESS serene rssSSeSs tsar SssSssorc rs 
‘ subroutine to calculate the fatigue stress concentration factor 
: based on Neuber’s method 
SUB Neuberkf 
CALL TYPIN("Enter notch sensitivity factor (q; 0 if unknown)", q) 
IF q <> 0 THEN é 
Ktf = 1+q* (Kt - 1) 
ELSE 
CALL TYPIN("Enter the notch root radius (r)", r) 
CALL TYPIN("Enter the appropriate material constant (rho)", rho) 
Kt£ = 1+ (Kt - 1) / (1 + (rho / xr) * .5) 
END IF 


END SUB 


R. SUBROUTINE NEWMAT 


' sessresenseeesesereeeeeeSe2=== NEWMAT SSS SSS SSS SSS SSS SSS SS SSS SS S552 5ES25=5== 
: subroutine to manually enter material constants and 
. enter new materials into the materials database 


BER ST SBS S SS ST SSS SS SSS SH SHS ST SSS SSS SSS SHS SSS SSS SSS SS SSS SSSSSeSsesSeqSeserzSere=x5== 


SUB NEWMAT (Ky) 


CONST ENTER = 13, ESCAPE = 27 

CONST DOWNARROW = 80, UPARROW = 72, LEFTARROW = 75, RIGHTARROW = 77 
CONST COL1 = 10, COL2 = 50, ROW = 7 

CONST COL3 = 7, COL4 = 62 

CONST Fields = 14 


’ conversion constanc from "ksi" to "MPa" 
CONST ab = 6.89476 


SHARED Menu AS Optionsl1, save$ 
DIM forceunit AS STRING, Fld AS INTEGER 


IF (Menu.material = 1) THEN 


forceunit = "ksi)" 
ELSE 

forceunit = "MPa)" 
END IF 


CALL MATMENU 


. Update field values and position based on keystrokes 
DO 

Put cursor on field 

LOCATE ROW + Fld, COL4 + 2 





: Get a key and strip null off if it’s an extended code 
DO 
Key$ = INKEYS 
LOOP WHILE Key$ = "" 


Ky = ASC(RIGHTS (Key$, 1)) 


SELECT CASE Ky 
CASE ENTER 
' Chech for all appropriate data parameters 
: before allowing return to main option menu 
IF b = 0 OR c = O OR Sigff = 0 OR epsff THEN 
LOCATE 23, 10: PRINT "Insufficient material parameters were entered" 
END IF 
CASE UPARROW, DOWNARROW 
‘ Adjust field location 
IF Ky = DOWNARROW THEN Inc = 1 ELSE Inc = -l 
Fld :- Rotated(0, Fields, Fld, Inc) 
CASE RIGHTARROW, LEFTARROW 
’ Adjust field value 
IF Ky = RIGHTARROW THEN Inc = 1 ELSE Ince = -1 
SELECT CASE Fld 


CASE 0 
‘ Ultimate Strength 
CALL TYPIN("Enter Ultimate Strength (Su)", Su) 
CALL MATMENU 


CASE 1 
' Yield Strength 
CALL TYPIN("Enter Yield Strength (Sy)", Sy) 


CALL MATMENU 


CASE 2 
‘ Fatigue Yield Strength 
CALL TYPIN("Enter Fatigue Yield Strength (Sy’)", Syf) 
CALL MATMENU 


CASE 3 
' Strength Coefficient 
CALL TYPIN("Enter Strength Coeficient (K)", K) 
CALL MATMENU 


CASE 4 
’ Cyclic Strength Coeficient 
CALL TYPIN("Enter Cyclic Strength Coeficient (K’)", K£) 
CALL MATMENU 


CASE 5 
‘ Cyclic Strength Coeficient 
CALL TYPIN("Strian Hardening Exponent (n)", n) 
CALL MATMENU 


CASE 6 
’ Cyclic Strian Hardening Exponent 
CALL TYPIN("Cyclic Strian Hardening Exponent (n’')", nf) 
CALL MATMENU 





CASE 7 
’ Ductility Coefficient (epsilon f) 


CALL TYPIN("Ductility Coefficient (epsilon f)", epsf) 


CALL MATMENU 


CASE 8 
' Fatigue Ductility Coefficient (epsilon f’) 


CALL TYPIN("Fatigue Ductility Coefficient (epsilon f’)", 


epsff) 
CALL MATMENU 
CASE 9 
‘ Strength Coefficient (sigma f) 
CALL TYPIN("Strength Coefficient (sigma f)", sigf) 


CALL MATMENU 


CASE 10 
' Fatigue Strength Coefficient (sigma f’) 


CALL TYPIN("Fatigue Strength Coefficient (sigma f’)", 


CALL MATMENU 


CASE 11 
' Patigue Strength Exponent (b) 
CALL TYPIN("Patigue Strength Exponent (b)", b) 
CALL MATMENU 


CASE 12 
' Fatigue Ductility Exponent (c) 
CALL TYPIN("Fatigue Ductility Exponent (c)", c) 
CALL MATMENU 


CASE 13 
' Endurance Strength (Sf) 
CALL TYPIN("Endurance Strength (Sf)", Sf) 
CALL MATMENU 


CASE 14 
' Modulus of Elasticity (E) 
CALL TYPIN("Modulus of Elasticity (E)"”, E&) 
CALL MATMENU 


CASE ELSE 
END SELECT 


CASE ELSE 
END SELECT 


‘ exit do loop and continue with subroutine if ENTER 
LOOP UNTIL (Ky = ENTER OR Ky = ESC) 


CLS 
LOCATE 13, 5: PRINT "Enter Material’s Name (up to 30 charactors): " 
LOCATE 13, 52: INPUT matname$ 


LOCATE 15, 5: PRINT "Save this material in material data base (Y/N): " 
LOCATE 15, 54: INPUT saveS$ 
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sigff) 


‘ progam segment to save new materials in the material data base 
IF save$ = "YY" OR save$ = "y" THEN 

‘OPEN "A:\MAT.DAT" FOR APPEND AS #3 

‘OPEN "B:\MAT.DAT" FOR APPEND AS #3 

OPEN "MAT.DAT" FOR APPEND AS #3 


matcount = matcount + 1 
IF fForceunit = "ksi)" THEN 
matindex = matcount 
WRITE #3, matindex, matname$, Su, Sy, Syf, K, Kf, n, nf, epsf, epsff, 
sigf, sigff, b, c, Sf, SfSu, E 
WRITE #3, matindex + 50, matname$, Su * ab, Sy * ab, Syf * ab, K * ab, 
Kf * ab, n, nf, epsf, epsff, sigf * ab, sigff * ab, b, c, Sf 
* ab, SfSu, E * ab 
ELSE 
matindex = matcount + 50 
WRITE #3, matindex, matname$, Su, Sy, Syf, K, Kf, n, nf, epsf, epsff, 
sigf, sigff, b, c, Sf, Sf£Su, E 
WRITE #3, matindex - 50, matnameS, Su / ab, Sy / ab, Syf / ab, K / ab, 
Kf / ab, n, nf, epsf, epsff, sigf / ab, sigff / ab, b, c, Sf 
/ ab, SfSu, E / ab 
END IF 


CLOSE #3 


‘OPEN "A:\NEWCOUNT.DAT" FOR OUTPUT AS #4 
‘OPEN "B:\NEWCOUNT.DAT" FOR OUTPUT AS #4 
OPEN "NEWCOUNT.DAT" FOR OUTPUT AS #4 
newstuff = newstuff + 1 

WRITE #4, newstuff 

CLOSE #4 


END IF 


CLS 
CALL UPDATEMENU 
Ky = 0 


END SUB 


S. SUBROUTINE OPTMENU1 


SSSSS2SSsseSSSsessqsqeqssaqqeo2zx Option Menu 1 BSSsssssessesesaesssescseessrsezec2z= 
r Define/swtich the user selectable functions for the material and 
: stress concentration factors 


SUB OPTMENU1 STATIC 


SHARED VC AS Config, Menu AS Optionsl, Available AS STRiNG, Fields AS INTEGER 
SHARED Fld AS INTEGER 


‘ Constants for key codes and column positions 

CONST ENTER = 13, ESCAPE = 27, Fl = 59 

CONST DOWNARROW = 80, UPARROW = 72, LEFTARROW = 75, RIGHTARROW = 77 
CONST COL1i = 10, COL2 = 50, ROW = 9 

CONST COL3 = 7, COL4 = 42 
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‘OPEN "A:\MAT.DAT” FOR INPUT AS #3 

‘OPEN "B:\MAT.DAT" FOR INPUT AS #3 

OPEN "MAT.DAT" FOR INPUT AS #3 

INPUT #3, matcount : 
CLOSE #3 

‘OPEN "A:\NEWCOUNT.DAT" FOR INPUT AS #4 

‘OPEN "B:\NEWCOUNT.DAT" FOR INPUT AS #4 . 
OPEN "NEWCOUNT.DAT" FOR INPUT AS #4 

INPUT #4, newstuff 

CLOSE #4 

matcount = matcount + newstuff 

CALL Loadmaterial (Menu.material) 


‘ Block cursor 
LOCATE ROW, COL1, 1, 1, 12 


CALL UPDATEMENU 


Skip field 10 if there’s only one value 
IF LEN(AvailableS) = 1 THEN Fields = 8 ELSE Fields = 10 


' Update field values and position based on keystrokes 
DO 

£ Put cursor on field 

LOCATE ROW + Fld, COL4 + 2 


f Get a key and strip null off if it’s an extended code 
DO 
Key$ = INKEYS 
LOOP WHILE Key$ = "" 


Ky = ASC(RIGHTS (Key$, 1)) 


SELECT CASE Ky 
CASE ESCAPE 
‘ End program 
CLS : END 
CASE F1 ’ review material parameters 


DO 

Key2$ = INKEYS 

LOOP WHILE Key2$ = "" 
Ky2 = ASC(RIGHTS (Key2$, 1)) 
LOOP UNTIL Ky2 = ESCAPE 


CLS : 
CALL UPDATEMENU 
CASE ENTER 
’ runs material parameter entry subroutine 
IF Menu.material = 1 OR Menu.material = 51 THEN NEWMAT (Ky) . 


CASE UPARROW, DOWNARROW 
’ Adjust field location 
IF Ky = DOWNARROW THEN Inc = 2 ELSE Inc = -2 
Fld = Rotated(0, Fields, Fld, Inc) 
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CASE RIGHTARROW, LEFTARROW 
‘ Adjust field value 
IF Ky = RIGHTARROW THEN Inc = 1 ELSE Inc = -l 
SELECT CASE Flid 


CASE 0 
, ‘ Units 
IF Menu.units = "Brit" THEN 
Menu.units = " SI " 


: Menu.material = Menu.material + 50 
CALL Loadmaterial (Menu.material) 
CALL UPDATEMENU 


Menu.units = "Brit" 

Menu.material = Menu.material - 50 
CALL Loadmaterial (Menu.material) 
CALL UPDATEMENU 


END IF 
CASE 2 
' Material Selection 
IF Menu.units = "Brit" THEN 


Menu.material = Rotated(1, matcount, Menu.material, Inc) 
CALL Loadmaterial (Menu.material) 

PRINT USING "&"; matname 

ELSE 
Menu.material = Rotated(51, matcount + 50, Menu.material, 
Inc) 
CALL Loadmaterial (Menu.material) 
PRINT USING "&"; matname 


END IF 
CASE 4 
‘ Stress Consentration Factor Kt 
CALL TYPIN("Enter stress concentration factor (Kt): ", Kt) 


CALL UPDATEMENU 


CASE 6 
‘ Patigue Stress Concentration Calculation Method 
Menu.Ktf = Rotated(i, 3, Menu.Ktf, Inc) 
SELECT CASE Menu.kKtf 


CASE 1 
' manual entry‘ 
option4 = "Manually Entered " 
LOCATE ROW + 8, COL4: PRINT USING "&"; "[ --- ] " 
LOCATE ROW + 6, COL4 + 2: PRINT USING "&"; option4 
CASE 2 
' Neubers Method 
option4 = "Neuber’s Method " 
LOCATE ROW + 8, COL4: PRINT USING "&"; "[ --- ] Bah 
LOCATE ROW + 6, COL4 + 2: PRINT USING "&"; option4 
CASE 3 
‘ Peterson Method 
, option4 = "Peterson’s Method" 
LOCATE ROW + 8, COL4: PRINT USING "&"; "[ --- ] st 
LOCATE ROW + 6, COL4 + 2: PRINT USING "&"; option4 
CASE ELSE 


END SELECT 








CASE 8 
' Fatigue Stress Concentration Factor displayed 
SELECT CASE Menu.Ktf 
CASE 1 
‘ manual entry' 
CALL TYPIN("Enter fatigue stress concentration factor 
{(Ktf): ", Ktf) 
CALL UPDATEMENU 
CASE 2 
'  Neubers Method 
CALL NeuberKf 
CALL UPDATEMENU 
CASE 3 
' Peterson Method 
CALL Petersonkf 
CALL UPDATEMENU 
CASE ELSE 
END SELECT 


CASE 10 
' Available screen modes 
i = INSTR(Available$, HEX$(VC.Scrmn) ) 
i = Rotated(1, LEN(Available$), i, Inc) 
vc.Sern = VAL("&h" + MIDS (Availables, i, 1)) 
PRINT USING "##"; vVC.Scrn 

CASE ELSE 

END SELECT 


CASE ELSE 
END SELECT 


' aif cluase to ensure selection of a fatigue stress concentation factor 
IF Ky = ENTER AND Ktf = 0 THEN 


Fld = 8 
CALL TYPIN("Enter fatigue stress concentration factor (Ktf): ", Ktf) 
CALL UPDATEMENU 
Ky = 0 
END IF 


‘ Return to main program if ENTER 
LOOP UNTIL (Ky = ENTER AND Ktf <> 0 AND matname <> "new material”) 


END SUB 
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T. SUBROUTINE OPTMENU2 


’ St Option Menu 2 SSeseSeSersertrerssssssssserrerereses 
, Define/swtich the user selectable functions for processing options 


SUB OPTMENU2 STATIC 


SHARED VC AS Config, menu2 AS Options2, Available AS STRING, Fields AS INTEGER 
SHARED Fld AS INTEGER 


‘ Constants for key codes and column positions 

CONST ENTER = 13, ESCAPE = 27, F2 = 60 

CONST DOWNARROW = 80, UPARROW = 72, LEFTARROW = 75, RIGHTARROW = 77 
CONST COL1l = 10, COL2 = 50, ROW = 9 

CONST COL3 = 7, COL4 = 42 


’ Return cursor to menu top 
Fld = 0 


CALL UPDATEMENU2 


4 Skip field 10 if there’s only one value 
IF LEN(Available$) = 1 THEN Fields = 8 ELSE Fields = 10 


4 Update field values and position based on keystrokes 
DO 

i Put cursor on field 

LOCATE ROW + Fld, COL4 + 2 


: Get a key and strip null off if it’s an extended code 
DO 

Key$ = INKEYS$ 

LOOP WHILE Key$ = "" 


Ky = ASC(RIGHTS(Key$, 1)) 
SELECT CASE Ky 


CASa ESCAPE 
‘ End program 


r CLS : END 
CLS : flag2 = 0 
CASE F2 ’ changes sound to condition show in flag5 
IF flagS = "ON" THEN ‘ (sound is opposite of flag) 
flags = "OFF" 
ELSE 
flagS = "ON" 
END IF 
UPDATEMENU2 


CASE UPARROW, DOWNARROW 
‘ Adjust field location 
IF Ky = DOWNARROW THEN Inc = 2 ELSE Inc = -2 
Fld = Rotated(0, Fields, Fld, Inc) 

CASE RIGHTARROW, LEFTARROW 
‘ Adjust field value 
IF Ky = RIGHTARROW THEN Inc = 1 ELSE Inc = -1l 
SELECT CASE Fld 
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CASE 0 
‘ type of inputs (stress vs strain) 
IF menu2.inputs = "Stress" THEN 
menu2.inputs = "Strain" 
inputfile = "STRAIN.DAT" 
ELSE 
menu2.inputs = "Stress" 
inputfile = "STRESS.DAT" 
END IF 
UPDATEMENU2 
CASE 2 


‘ Equation Selection 
menu2.equation = Rotated(1, 3, menu2.equation, 
SELECT CASE menu2.equation 


CASE 1 
' Morrow's equation 
equationname = "Morrow's equation" 
CASE 2 
' Smith Watson Topper 
equationname = "Smith-Watson-Topper" 
CASE 3 
‘ Manson Halford 
equationname = "Manson-Halford" 
CASE ELSE 
END SELECT 
UPDATEMENU2 


CASE 4 
' variable/fixed n’ and K’ 


Inc) 


menu2.option3 = Rotated(1, 4, menu2.option3, Inc) 


SELECT CASE menu2.option3 


CASE 1 

’ fixed n and K 

nKtype = " Fixed n’ and K’ " 
CASE 2 


‘ variable n’ and XK’ 
nKtype = "Variable n’ and K'" 
CASE 3 
' variable n’ and fixed K’ 
nKtype = "Variable n’ and fixed K’" 
CASE 4 
' fixed n’ and variable XK’ 
nKtype = "Fixed n’ and variable K’'" 
CASE 5 
' Experimental - not used 
nKtype = "Experimental" 
CASE ELSE 
END SELECT 
UPDATEMENU2 


CASE 6 
‘ calculation type 


menu2.option4 = Rotated(1, 3, menu2.option4, Inc) 


SELECT CASE menu2.option4 


CASE 1 

‘ blocks to failure 

calctype = "Load blocks to failure" 
CASE 2 

' single block effects 

calctype = "Single block effects" 
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CASE 3 
' batch process 
Ccalctype = "Batch Process" 
flag4 = 1 
CASE ELSE 
END SELECT 
UPDATEMENU2 


CASE 8 
‘ input data file name 
CALL TYPINSTRING ("Input data file’s name:", 


UPDATEMENU2 


inputfile) 


CASE 10 
‘output data file name 
CALL TYPINSTRING ("Output data file’s name:", outputfile) 


UPDATEMENU2 


CASE ELSE 
END SELECT 


CASE ELSE 
END SELECT 


‘ if cluase to ensure selection of a fatigue stress concentation factor 


IF Ky = ENTER AND Ktf = 0 THEN 


Fld = 8 
CALL TYPIN("Enter fatigue stress concentration factor (Ktf): ", Ktf) 


CALL UPDATEMENU 
Ky = 0 
END IF 


* Return to main program if ENTER 


LOOP UNTIL (Ky = ENTER AND Ktf <> 0 AND matname <> "new material") OR (Ky = 
ESCAPE) 


CLS 


END SUB 
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U. SUBROUTINE OUTPUTER 


eS ee ee et ee ee ee ee ee ee et OUTPUTER eee eee ee ee er ee ee ee Pe ee ee 
: subroutine writes to an output file compution parameters 
and results 


SUB OUTPUTER 
OPEN outputfile FOR APPEND AS #11 


WRITE #11, 

WRITE #11, 

WRITE #11, equationname 

WRITE #11, nKtype 

WRITE #11, calctype 

WRITE #11, “input file:", inputfile 
WRITE #11, "output file:", outputfile 
WRITE #11, "blocks:", block 

WRITE #11, "i counter:", lasti 

WRITE #11, "reversal count:", NNfcount 
WRITE #11, "life factor:", usedlife 
WRITE #11, 

WRITE #11, 


CLOSE #11 


END SUB 


Vv. SUBROUTINE PetersonKf 


SSsese22se222 SS 2525525555555 Petersonk£ SSSsese2eSSSESeSSseeeseqorsszsses 
subroutine to calculate the fatigue stress concentration factor 
based on Peterson’s method 


SUB PetersonKf 


CALL TYPIN("Enter notch sensitivity factor (q; 0 if unknown)", q) 
IF q <> 0 THEN 
Ktf = 1+q* (Kt - 1) 
ELSE 
CALL TYPIN("Enter the notch root radius (r)", r) 
CALL TYPIN("Enter the appropriate material constant (a)", a) 
Ktf£ = 1 + (Kt - 1) / (1 +a/-r) 
END IF 
END SUB 
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W. FUNCTION Rotated 


: tt? 2228 _2_8_-$_-8_-$--S—-4-2-4--5_-3 5-5 Rotated Serres sesseseeseesrer res TreSeSeSsSseSaeese=== 
: Returns the present value adjusted by Inc and rotated if necessary 
. so that it falls within the range of Lower and Upper. 


' Calculate the next value 
present = present + Inc 


' Handle special cases of rotating off top or bottom 
IF present > Upper THEN present = Lower 

IF present < Lower THEN present = Upper 

Rotated = present 


END FUNCTION 


X. SUBROUTINE TYPIN 


j subroutine to update the value of a specific variable through 
keyboard entry (real number variable) 

SUB TYPIN (cstring$, valu) 

LOCATE 23, 10: PRINT cstring$ 

LOCATE 23, 58: INPUT valu 

LOCATE 23, 10: PRINT " 

" 

END SUB 

Y. SUBROUTINE TYPINSTRING 

§ sere re ress esseseeeeeees=ea== TYPINSTRING pmseestet esses SeeeeeeeeeSeese= 


subroutine to update the value of a specific string variable through 
’ keyboard entry 


SUB TYPINSTRING (cstring$, stringname$) 
LOCATE 23, 10: PRINT cstring$ 

LOCATE 23, 58: INPUT stringname$ 
LOCATE 23, 10: PRINT " 


END SUB 
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Z. SUBROUTINE UPDATEMENU 


' subroutine to update the options menu (menu 1) 


SUB UPDATEMENU 


SHARED VC AS Config, Menu AS Optionsl, Available AS STRING, 


Fld AS INTEGER 


‘ Constants for key codes and column positions 
CONST COL1 = 10, COL2 = 50, ROW = 9 
CONST COL3 = 7, COL4 = 42 


’ Display key instructions 
LOCATE 1, COL1 


PRINT "OB. seca cw uetes,’s oad Move to next field" 
LOCATE 2, COL1 
PRINT "DOWN ........ Move to previous field" 
LOCATE 3, COLL 
PRINT "LEFT/RIGHT .... Change field up/down" 
LOCATE 4, COL1 
PRINT "Fl .... Display matrial’s parameters" 
LOCATE 5, COL1 
PRINT "ENTER .... Start with current values” 
LOCATE 6, COL1 
PRINT "ESCAPE ................ Quit Program" 


’ Display fields 


LOCATE ROW, COL3: PRINT "Type of units (SI or British)"; 
LOCATE ROW, COL4: PRINT USING "[{ & ]"; Menu.units; 


LOCATE ROW + 2, COL3: PRINT "Material"; 


Fields AS INTEGER, 


LOCATE 


LOCATE 
LOCATE 


LOCATE 
LOCATE 


LOCATE 
LOCATE 


LOCATE 
LOCATE 


END SUB 


ROW + 2, COL4: 


ROW 
ROW 


ROW 
ROW 


ROW 
ROW 


ROW 
ROW 


+ 
+ 


+ 
+ 


+ 


4, COL3: 
4, COL4: 
6, COL3: 
6, COL4: 
8, COL3: 
8, COL4: 
10, COL3: 
10, COL4: 


PRINT 


PRINT 
PRINT 


PRINT 
PRINT 


PRINT 
PRINT 


USING "[ & ]"; matname; 


"Stress concentration factor (Kt)"; 
USING "[( ##.### ]"; Kt; 


"Method to calculate Kf"; 
USING "{ & ]"; option4 


"Fatigue stress conc. factor (Kf)"; 
USING "[( ##.### ]"; Ktf; 


PRINT "Screen Mode"; 
PRINT USING "[ ## ]"; VC.Scern 
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AA. SUBROUTINE UPDATEMENU2 


§ sesescscSssesseSSee5eSeeeee=== UPDATEMENU2 St StS 


: subroutine to update the options menu (menu 2) 


SUB UPDATEMENU2 


SHARED VC AS Config, menu2 AS Options2, Available AS STRING, 
Fld AS INTEGER 


‘ Constants for key codes and column positions 


CONST COL1 = 10, COL2 = 50, ROW = 9 
CONST COL3 = 7, COL4 = 42 
CLS 


’ Display key instructions 
LOCATE 1, COL1 


PRINT: "UP 2. .c. eee sb eae s Move to next field" 
LOCATE 2, COL1 

PRINT "DOWN ........ Move to previous field" 
LOCATE 3, COL1 

PRINT "LEFT/RIGHT .... Change field up/down" 
LOCATE 4, COL1 

PRINT USING "F2 ............-.. Turn sound to & "; flag5 
LOCATE 5, COL1 

PRINT "ENTER .... Start with current values" 
LOCATE 6, COL1 

PRINT "ESCAPE ... Return to previous screen" 


' Display fields 


Fields AS INTEGER, 


LOCATE ROW, COL3: PRINT "Type of inputs (stress or strain)"; 


LOCATE ROW, COL4: PRINT USING "[ & ]"; menu2.inputs; 


LOCATE ROW + 2, COL3: PRINT "Equation"; 
LOCATE ROW + 2, COL4: PRINT USING "[({ & J"; equationname; 


LOCATE ROW + 4, COL3: PRINT "Fixed / Varing n’ and K’"; 
LOCATE ROW + 4, COL4: PRINT USING "[ & J"; nKtype; 


~ s 


LOCATE ROW + 6, COL3: PRINT "Calculation Type oo 
LOCATE ROW + 6, COL4: PRINT USING "{ & ]"; calctype; 


ss 


LOCATE ROW + 8, COL3: PRINT "Input file name uee 


LOCATE ROW + 8, COL4: PRINT USING "[ & ]"; inputfile; 
LOCATE ROW + 10, COL3: PRINT "Output file name "; 
LOCATE ROW + 10, COL4: PRINT USING "[ & ]"; outputfile; 
END SUB 
67 


AB. FUNCTION xxKf 


4 2 — + $44 — 4-55 —4 5-5 5 5-5 2 a xxKf — 2 
Returns the value calculated for K’ based on the 
u number of cycles executed ‘ 
FUNCTION xxKf 
SHARED menu2 AS Options2 é 
SELECT CASE menu2.option3 
CASE 1 'fixed K’ 
xxKf = Kf 
CASE 2 ' variable K’ 
IF NNfcount > 1000000 THEN 
xxKf = KE 
ELSEIF NNfcount < 2000 THEN 
xxKf = K 
ELSE 
xxKf = ((Kf - K) / (LOG10(500000) - LOG10(1000))) * 
(LOG10 (NNfcount / 2) - LOG10(1000)) + K 
END IF 
CASE 3 ‘fixed K’ 
xxK£f = Kf 
CASE 4 ’ variable K’ 
IF NNfcount > 1000000 THEN 
xxKf = Kf 
ELSEIF NNfcount < 2000 THEN 
xxKf = K 
ELSE 


xxKf = ((K£ - K) / (LOG10(500000) - LOG10(1000))) * 
(LOG10(NNfcount / 2) - LOG10(1000)) + K 
END IF 
CASE ELSE 
END SELECT 


END FUNCTION 


AC. FUNCTION xxnuf 


7 > —— 2 — 225-4 4-5 4-5 5 4-3-5 4-2-5 2-5-5 5-5-1 + 5-1-4 4-3 xxnf I 
t Returns the value calculated for n’ based on the 

: number of cycles executed 

FUNCTION xxnf 

SHARED menu2 AS Options2 


SELECT CASE menu2.option3 


CASE 1 ‘ fixed n’ 
xxnf = nf 
CASE 2 ‘ variable n’ 


IF NNfcount > 1000000 THEN 
xxnf = nf 








ELSEIF NNfcount < 2000 THEN 


xxnf =n 
ELSE 
xxnf = ((nf - n) / (LOG10(500000) - LOG10(1000))) * 
(LOG10 (NNfcount / 2) - LOG10(1000)) +n 
END IF 
CASE 3 ’ variable n’ 


IF NNfcount > 1000000 THEN 
xxnf = nf 
ELSEIF NNfcount < 2000 THEN 


xxnf =n 
ELSE 
xxnf = ((nf - n) / (LOG10(500000) - LOG10(1000))) * 
(LOG10 (NNfcount / 2) - LOG10(1000)) +n 
END IF 
CASE 4 ' fixed n’ 
xxnf = nf 
CASE ELSE 


END SELECT 


END FUNCTION 


AD. FUNCTION xxNNfcount 


4 seeccsescsssesseesssesesee= XXNNfcount =<===sss=sssssssssess=s5e5snee55 
f Returns the number of cycles executed 


SST SSSSSSSSSSS SSS SST SSS SH SS SSS HSS SH SS SSS SSS SSS SS SS SSS SH SH SS SS HSS SSS SST T= 


FUNCTION xxNNfcount& 
‘DIM xxNNfcount AS LONG 


xxNNfcount = (top * block + i) 
END FUNCTION 
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APPENDIX B. PROGRAM LOADGEN 


Random load generation program 


es ss 


based on a typicical 1000 hour block for an A-6 aircraft 


DIM outfile AS STRING 
CONST fourg 1978 
CONST fiveg 333 
CONST sixg = 48 
CONST seveng = 10 
CONST totalg = fourg + fiveg + sixg + seveng 
CONST Su = 84 

CONST Sy = 76 

CONST gdesign = 6.5 


RANDOMIZE TIMER 


CLS 
‘LOCATE 20, 10: PRINT "Enter name for stress load output file:" 
‘LOCATE 20, 50: INPUT outfile 


CLS 

FOR j = 1T04 *’ loop to create four random files 
fourcount = fourg 

fivecount = fiveg 

sixcount = sixg 

sevencount seveng 


totalcount = fourg + fiveg + sixg + seveng 


SELECT CASE j 


CASE 1 

OPEN "testaa" FOR OUTPUT AS #2 
CASE 2 

OPEN "testbb" FOR OUTPUT AS #2 
CASE 3 

OPEN "testcc" FOR OUTPUT AS #2 
CASE 4 

OPEN "testdd" FOR OUTPUT AS #2 
CASE ELSE 


END SELECT 





"g" load history greneration: 
OPEN "gseries" FOR OUTPUT AS #1 


, WHILE totalcount > 0 
DO 
x = RND 


IF x <= (seveng / totalg) THEN 
7 7+ "g" case 
IF sevencount = 0 THEN EXIT DO 
y = 7 + (INT(xx * 10)) / 10 
WRITE #1, y 
WRITE #1, 1 
sevencount = sevencount - 1 
ELSEIF x <= ((seveng + sixg) / totalg) THEN 
‘ 6 to 7 "g" case 
IF sixcount = 0 THEN EXIT BO 
y = 6 + (INT(xx * 10)) / 10 
WRITE #1, y 
WRITE #1, 1 
sixcount = sixcount - 1 
ELSEIF x <= ((seveng + sixg + fiveg) / totalg) THEN 
’ 5 to 6 "g" case 
IF fivecount = 0 THEN EXIT DO 
y= 5 + (INT(xx * 10)) / 10 


WRITE #1, y 

WRITE #1, 1 

fivecount = fivecount - 1 
ELSE 


' 4 to 5 "g" case 

IF fourcount = 0 THEN EXIT DO 

y = 4 + (INT(xx * 10)) / 10 

WRITE #1, y 

WRITE #1, 1 

fourcount = 
END IF 


fourcount - 1 


totalcount = fourcount + fivecount + sixcount + sevencount 


LOOP 
WEND 
WRITE #1, 999 
CLOSE #1 





conversion of "g" load history to a stress load history: 


IF Su > (1.5 * Sy) THEN 
gtostress = Su / (1.5 * gdesign) 
ELSE 
gtostress = Sy / (qdesign) 
END IF 


OPEN "gseries"” FOR INPUT AS #3 
INPUT #3, load 
WHILE load < 9399 

stressload = load * gtostress 


WRITE #2, stressload 
INPUT #3, load 


WEND 
CLOSE #2 
CLOSE #3 
load = 0 
NEXT j 
END 
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APPENDIX C. MATLAB DATA REDUCTION PROGRAM 


A. PROGRAM CODE 


% UNCORRECTED STRAIN + LOAD DATA (for K and n) 

format compact 

% Monotonic #1 

% (81 points) 

% 0 1 2 3 4 5 6 7 8 9 

$ sera seeSS SESH SSS HS SSH SS SH SSS SS SSS SSS SS SS SS SS SS SS SS SS SSS SSS ese SSS Sess Ssscss 


ldOa=[ 305; 628; 998; 1332; 1650; 1979; 2307; 2627; %2975; 3350; 


17884; 17921; 17967; 17990; 18027; 18060; 18266; 18419; 18532; 18627; 
18901; 18843; 18916; 18932; 19065; 19128; 19176; 19225; 19264; 19302; 


% “ionotonic #2 
% (88 points) 


WwW 
ie] 
ul 
@ 
e 
my 
uy 
wo 
as 
WwW 
ul 
Oo 
> 
vi 
roa) 
=) 
nm 
~] 
Nn 
nN 
» 
\o 
fon) 
un 
wi 
b 
an 
an 
uw 
w 
ten) 
foe) 
uy 
wi 
~ 
NO 
paar 


ldOb=[ 306; 619; 962; 1303; 1619; 1958; 2275; 2590; 2911; 3242; 


18851; 18927; 19002; 19052; 19105; 19153; 19200; 19241; 19277; 19305; 
19224; 19258; 19381; 19404; 19417; 19433; 19438; 19449; 19452]; 
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% 10 percent: 3790 cycles 
% (73 points) 


% 0 1 2 3 4 5 6 T 8 9 

SSS 5SSSSs5SSsS255 SS SS 2S>S22535 52255 >S25S=SsSS222=>35 S555 5S2=S25>SSS 2255552555522 5522S S2= 

ex10=[.0005; .0013; 0022; .0030; .0038; .0047; .0054; .0063; .0071; .007¢ 
.0088; 0099; 0208; 0117; 0126; 0134; 0142; 0151; 0159; 0167 
-0175; .0184; .0192; .0200; .0209; .0217; .0225; .0233; .0242; .0250 
.0258; .0266; .0274; .0283; .0291; .0299; .0307; .0316; .0324; .0332 
.0361; .0381; .0402; .0422; .0444; .0466; .0489; .0511; .0532; .0553 
.0574; .0595; .0618; .0638; .0660; .0682; .0703; .0725; .0746; .0778 
.0798; .0818; .0840; .0861; .0882; .0902; .1102; .1304; .1506; .1707 
-1911; .2114; 2314; .2516]; 

1dio0=[ 302; 595; 934; 1252; 1561; 1864; 2160; 2469; 2770; 3098 
3446; 3795; 4170; 4486; 4785; 5076; 5365; 5671; 5967; 6274 
6571; 6856; 7161; 7450; 7748; 8026; 8334; 8619; 8911; 9211 
9490; 9787; 10071; 10361; 10652; 10946; 11225; 11517; 11812; 12074 
13059; 13751; 14478; 15148; 15875; 16591; 17315; 18004; 18577; 19073 
19372; 19536; 19560; 19573; 19577; 19572; 19553; 19524; 19483; 19365 
19325; 19297; 19303; 19297; 19282; 19324; 19304; 19264; 19261; 19243 
19202; 19137; 19015; 18802]; 

% 20 percent: 7580 cycles 

$ (87 points) 

% G 1 2 3 4 5 6 7 8 9 

ee ee ee ee ee 

ex20=[.0007; .0016; .0028; .0039; .0049; .0060; .0069; .0081; .0090; .0162 
.0113; .0121; .0132; .0142; .0153; .0162; .0172; .0181; .0191; .02U1 
.0211; .0220; .0229; .0240; .0250; .0260; .0271; .0282; .0295; .0309 
.0322; .0334; .0349; .0363; .0379; .0394; .0410; .0426; .0442; .0458 
.0515; .0536; .0558; .0581; .0604; .0626; .0649; 0672; .0695; .0717 
.0740; .0762; .0784; .0806; .0826; .0849; .0873; .0897; .0918; .1010 
.1031; .1051; 1073; 1094; Ll. 7; 1139; 1339; 1542; 1747; 1948 
-2150; .2351; .2553; .2756; .2959; .3159; .3360; .3561; .3764; .3965 
4168; .4368; .4569; .4771; .4974; .5179; .5385; .5586); 

ldz0=[ 327; 650; 984; 1324; 1641; 1955; 2271 2609; 2952; 3291 
3592; 3911; 4246; 4540; 4896; S206; 5513; 5827; 6132; 6468 
6782; 7077; 7411; 7699; 8034; 8301; 8629; 8918; 9245; 9516 
9816; 10093; 10364; 10670; 10939; 11241; 11518; 11793; 12051; 12331 
13238; 13543; 13876; 14198; 14513; 14846; 15122; 15467; 15741; 16030 
16302; 16577; 16838; 17029; 17194; 17366; 17483; 17593; 17656; 17921 
17944; 18007; 18062; 18086; 18120; 18157; 18402; 18583; 18703; 18802 
18888; 18951; 19025; 19092; 19140; 19209; 19253, 19297; 19343; 19387 
19414; 19456; 19463; 19501; 13528; 19576; 19567; 19571); 





AP AP 


° 
Qa 


ex30=[.0004; 
.0089; 
.0174; 
.0255; 
.0355; 
.0568; 
.0843; 
.1989; 
.4016; 


1d30=[ 293; 
3546; 
6812; 
9809; 
13473; 
19101; 
19093; 
19181; 
19489; 


AP AP 


30 percent: 
(88 points) 


40 percent: 
(86 points) 


19073; 
12239; 
19498; 


19093; 
19114; 
19255; 
19501; 


11370 cycles 


10719; 
15811; 
19094; 
19104; 
19295; 
19493; 


15160 cycles 


19438; 


19498]; 


.1780; 
.3813; 


3196; 
6486; 
9520; 
12473; 
19085; 
19089; 
19132; 
19486; 


ae AO 
fo) 


ex40=[(.0002; 
.0129; 
.0286; 
.0460; 
.0697; 
.0911; 
.1214; 
.2343; 
.4369; 


1d40=[ 158; 
3194; 
6077; 
8836; 
9823; 
12099; 
17423; 
18920; 
19503; 


19008; 
19539; 


14430; 


17247; 
19074; 
19445; 
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19185; 
19664; 


19232; 
19667] ; 


18422; 
19362; 


1861° ; 
19395 


18789; 
19464; 


% adjust number: 


area = .046875; 

ld0Qa = ld0a / 5; stressOa = ld0a / area; 
1d0b = 1d0b / 5; stressOb = 1d0b / area; 
1d10 = 1d10 / 5; stress10 = 1d10 / area; 
1d20 = 1d20 / 5; stress20 = 1d20 / area; 
1d30 = 1d30 / 5; stress30 = 1030 / area; 
1d40 = 1d40 / 5; stress40 = 1d40 / area; 
strain0a = ex0a * .3; 

strainOb = ex0b * .3; 

strainlO = exl0 * .3; 

strain20 = ex20 * .3; 

strain30 = ex30 * .3; 

strain40 = ex40 * .3; 


$idel a:\ernie.met 
s!del ermie.met 


plot (strain0a, stress0a),title(’Monotonic #1 Stress-Strain’),grid 
xlabel(’Strain (in/in)’),ylabel(’Stress (psi)’); 

‘pause 

meta a:\ernie 


plot (strain0b, stress0b) ,title(’Monotonic #2 Stress-Strain’) ,grid 
xlabel(’Strain (in/in)'),ylabel(’Stress (psi)’); 

‘pause 

smeta a:\ernie 


plot (strainl0,stress10) ,title(’10% Life Stress-Strain’),grid 
xlabel(’Strain (in/in)’),ylabel(’Stress (psi)’); 

spause 

meta a:\ernie 


plot (strain20,stress20) ,title(’20% Life Stress-Strain’),grid 
xlabel(’Strain (in/in)’),ylabel(’Stress (psi)’); 

spause 

‘meta a:\ernie 


plot (strain30,stress30),title(’30% Life Stress-Strain’),grid 
xlabel(’Strain (in/in)’),ylabel('Stress (psi)’)j; 

‘pause 

meta a:\ernie 


plot (strain40,stress40) ,title(’40% Life Stress-Strain’) ,grid 
xlabel(’Strain (in/in)’),ylabel(’Stress (psi)'); 

tpause 

‘meta a:\ernie 


plot (strainOb, stress0b, strainl0,stress10,strain20,stress20,... 
strain30,stress30,strain40,stress40),title(’Stress-Strain’),grid 
xlabel('Strain (in/in)’),ylabel(’Stress (psi)’); 

‘pause 

meta a:\ernie 
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!del ernie.out 
diary erine.out 


% 


compute true strain from engineering 
strainQa) ; 
strainob) ; 
strainl0) ; 
strain20) ; 
strain30); 
strain40); 


tstrainoda 
tstrain0db 
tstrainlo 
tstrain20 
tstrain30 
tstrain40 


log(1 
log(1 
log(1 
log({1 
log(1 
log(1 


++ ett + 


2 


% compute true stress from engineering 
tstrain0a) ; 
tstrainob) ; 
tstrainl0) ; 
tstrain20) ; 
tstrain30) ; 
tstrain40) ; 


(1 
(1 
(1 
(1 
(1 
(1 


stress0a 
stress0Ob 
stress10 
stress20 
stress30 
stress40 


sigl 
sig2 
sig3 
sig4 
sig5S 
sig6 


++ eee 


* 
* 
* 
. * 
* 
* 


lsigl 
lsig2 
lsig3 
lsig4 
lsig5 
lsig6é 


1og10(sigl) ; 
10g10 (sig2) ; 
1logl10 (sig3) ; 
1o0g10 (sig4) ; 
1og10(sigS) ; 
1o0g10 {sig6) ; 


'oui wouet 


E 


10300000; 


x1 logl0(tstrain0a(2:82) - sigl(2:82)/E); 


x2 
x3 
x4 
xS 
x6 


lsigla 


1lsig5 
lsig6 


xxl 
xx2 
xx3 
xx4 
xx5 
xx6 


i 


logl10(tstrainOb - 
log10(tstrainl10o - 
log10(tstrain20 - 
logl10 (tstrain30 (2: 
1og10 (tstrain40 (2: 


lsigl (2:82); 
1sigS (2:89) ; 
lsig6 (2:87) ; 


a 
a 


polyfit (x1, 
polyfit (x2, 
polyfit (x3, 
polyfit (x4, 
polyfit (x5, 
polyfit (x6, 


lsig2 
lsig3 
lsig4 


lsigla, 


lsig5a, 
lsigé6a, 


sig2/E); 
sig3/E) ; 
sig4/E) ; 
89) - sig5(2:89)/E); 
87) - sig6(2:87)/E) ; 


1) 
1) 
1) 
1) 

1) 

1) 


’ 
’ 


’ 
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PROGRAM OUTPUT/RESULTS 


Kl = 10*(xx1(2)) 

K2 = 10% (xx2(2)) 

K3 = 10% (xx3(2)) 

K4 = 10*(xx4(2)) 

KS = 10% (xx5 (2) ) 

K6 = 10% (xx6(2)) 

nl = xxl1(1) 

n2 = xx2(1) 

n3 = xx3 (1) 

n4 = xx4(1) 

nS = xx5(1) 

n6é = xx6(1) 

diary 

B. 

xxX1 = 0.4425 
xXX2 = 0.4331 
xx3 = 0.6246 
xXx4 = 0.5063 
xx5 = 0.4261 
XX6 = 0.5902 
K1 = 4.0798e+005 
K2 = 3.8107e+005 
K3 = 1.1229e+006 
K4 = 4.3900e+005 
K5 = 3.6831e+005 
K6 = 4.8105e+005 
ni = 0.4425 

n2 = 0.4331 

n3 = 0.6246 

n4 = 0.5063 

nS = 0.4261 

n6 = 0.5902 


UT UT UT OV UT UT 


-6106 
-5810 
-0503 
-6425 
-5662 
-6822 
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APPENDIX D. 


A. BRITISH/AMERICAN UNITS 








MATERIAL DATA BASE 


Monotonic and Cyclic Strain Properties of Selected Engineering Alioys: American/British units 


S, 
Process Be. SuiSe Kik’ owey 2A = 10 
Matenal Descnption (ksi) tksivks) (ksvksi} nin en tksivksid b t (ksi) SS, 
Steel 
10S. 1009 HR. sheet xa IRA T2AT 0.16:0.32 P60 $2302 -0.1909 - 20 71 OSs 
1005S. {009° CLD. sheer fl) S836 76/7) O.N490 1) 1020.11 122/78 -0.07%3 -04] 2k 047 
1920 HLR. sheet 64 3R38 jo77 12 0.19018 0.9604) 103130  -O.12  -0.81 oa 034 
0OU0 Cast steel i es 34/46 0.30/0.12  0.62/0.28 109/95 -().082 -0.8] 2k 0.38 
Man-Ten H.R. sheet 74 STIS4 —14 0.20011 1.020.86 Ft8/117 -0.071 -0.68 38 O.St 
140 As forged 90 $86 _—— 0.22/0.18 0.93/0.61 182/223 -0.14 -0.87 as 0.28 
RQC-10) H.R. sheet 138 12R/R7 170208 0.06/0.14 1020.66 1937580 -0.07 -0.69 s9 G43 
4142 Drawn at temp 184 182/108 ~— —'0.18 0.38/0.22 1627210 -O.10 -0.51 44 a.28 
4142 Q&T 205 200'T20 —— 0.051/0.17 0.660 48 265265 -0.08 -(0.78 73 (36 
4142 Q&T 280 250'198 —— 0.048/0.13 0.43/0.09 318315 -0.08f -G6l RS ou 
4340) H.R. and annealed = 120 92/66 —— —!0.18 0.57/0.48 [88174 -0.098 -0.54 40 03 
4340 QaT yRO 170110 299. 0.066/0 14 0.84/0.73 240/240 -G.076 -0 62 Tl 040 
4340 Q&T 213 199:]20 -_-— —WAS 0.48/0.48 = 226/290 - 0.091 -a.60 68 0.32 
9262 Annealed $34 h6'76 783200 O27OtS Olote ISVISt -OO7b -0 47 Si 38 
9262 Q&T 148 tf4/9s —~197 igo? 44a) W777  =-P.O72 -0 60 &s 0 3K 
Aluminum 
1100.0 As teceived {6 149 —_—— —/O.TS8 2.0918 IR -0 706-069 . 033 
2024-T3 — 68 S862 A6/98 11.032/0.068 0.28/0.22 RVIEO -O.F24 -0.89 Pes 032 
2024-14 a 69 44/65 7 O.200 08 0420.21 9247-0 Td 082 as av 
5456-H3 _ 58 44'S? _—— 16 0 42/046 T4108 -O1l -067 18 Oat 
7078.16 _ R4 48°76 1206— 1.10146 0.41019 1OR/191 -O.426 -0.82 Ee 030 
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B. SI UNITS 


Monotonic and Cyclic Strain Properties of Selected Engineering Alloys: S! Units 








Ye 
Process Sy Lee ae ALK TET, HN = 107, 
Matenal Descnpuon (MPaa (MPa MPar «MPa MPar Ph en€, (MPa MPa h ‘ oMPat SS, 
Steel 
1005-1009 H.R. sheet 338 262228 E27 462 16012 ft oO 1  &sRadl “9 109-1 39 148 Ay 
1005-1009 C.D. sheet 14 400,248 S24290  VO4WO TT FOO PE  S4b S38 “0073 -a4l 9s Ya? 
1020 HILR. sheet 441 262241 “38772 OAID AR 096041 T1ORI6 O42 -ast 1§2 34 
0030° Cast steel 196 WAIT —— 30/002 OA2O28R Sas} -O 082-8) 190 38 
Man-Ten HER. sheet SEO IDR/AT2 — TRA 0.200.112 tL O2Z086 814807 -07E - 48 762 VS} 
1040 As forged 62) WAL IRA —— W.220.18 O9VO6F 1OSH1S40 -0 Fs - SF 173 38 
RQC-100 HR. sheet QP 8RVEOD LTT VIIA 060.14 1.020.66  1330/1290 0-007 -4Y 403 43 
41420 Drawn at temp 1062 HOSR/745 a — O18 OF8022 TSt4so -H 10 -H SY 310 28 
4142 Q&T 1413) 1379-827 —— O.050097 C660 45 1828/1828 -N0OR -a7s 603 0 3A 
4142 Q&T 1931 1724/9344 ae 9.048/0. 13) 0.43/0.09 2170/2170 -.08T -9 61 SR9 ut 
4340 H.R. and annealed = 827 $34,458 —_— — OR 0 $70.45 1090/1200 -0.098 -0 S54 4 Os 
4340 Q&T 4d TE72-78SR_— STH. O.0660.14 0840.73 165S8/tHSS -0 07h -0 42 49? 0 40 
4340 Q&T 1469 137DR27 cae — O18 1.480.48 18602000 -0091 -160 th7 32 
9262 Annealed 924 AS88594 1THL1379 0 I2OAS = OT6O 16 10461046 -077 -0 47 24k IR 
9262 Q&T 1000 TRAHAR ISR OPO SVN ST 12209220 -077 -0 60 381 IR 
Aluminum 
1100-0 AS received 1h 9762 dated —O 1S 2991.8 —.19% 106 069 17 WA 
2024-T3 ~ 469 -3:79/427 488688 0.03270. 065 0 28.0.22  SSR/TIOO -O 12s -0 59 {St 0 32 
2024-T4 ots 476 303/441 807" 20008 O4SF0.21  A3atOls -a Th -as? 17s it 
$456- 113 ~ 406 234/359 —— —U.16 O.420.46 $24,725 mit -ae7 124 | 
7075-T6 _— S79 AAG S24 R27 O10 146 04119 TH8IRTS  -0 126-082 176 90 
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